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ABSTRACT 

This  study  was  conducted  to  examine  the  efTects  of  discrete  precipitation  events  on 
the  short-term  and  seasonal  evolution  of  ocean  mixed  layer  temperature  and  salinity 
structure.  This  study  was  located  at  Ocean  Station  "P"  (50°N,  145°W)  in  the  Northeast 
Pacific  Ocean.  Two  numerical  modeling  experiments  were  performed.  The  first  was  to 
simulate  the  response  of  the  mixed  layer  to  hypothetical  discrete  (isolated)  precipitation 
events.  This  experiment  showed  that  the  effect  of  a  single  discrete  rain  event  can  van' 
with  season,  with  the  mixed  layer  depth  (MLD)  at  onset  of  the  rain  event,  and  with  the 
strength  of  wind  stress  forcing.  A  single  rain  event  can  have  lasting  effects  on  mixed 
layer  depth  and  temperature  for  up  to  55  days  after  the  event,  depending  upon  the  sea- 
son. 

The  second  experiment  simulated  quasi-realistic  "complex"  precipitation  forcing, 
with  a  realistic  distribution  of  synoptic  events  over  a  13-month  period.  For  this  exper- 
iment, four  different  precipitation  "intensities"  were  tested.  The  values  of  temperature, 
salinity,  and  MLD  predicted  by  the  model  were  compared  with  observed  Conductivity 
Temperature  Depth  measurements  and  with  the  values  predicted  using  constant  precip- 
itation forcing.  In  all  experiments,  the  modelled  MLD's  approximated  the  observed 
MLD  and  temperature  cycle.  MLD's  for  all  cases  were  too  deep  while  temperature  was 
estimated  well  in  winter  and  was  too  cool  in  summer.  Predicted  salinity  provided  the 
greatest  discrepancy  between  the  modeled  and  the  observed  cycle.  The  96  cm,  400-day 
quasi-realistic  precipitation  forced  case  best  approximated  the  salinity  observations 
though  predicted  salinity  was  fresher  than  observed  in  winter  and  saltier  than  observed 
in  summer.  Model  results  suggest  that  the  amount  of  precipitation  observed  at  Ocean 
Station  "P"  was  too  low  to  explain  the  observed  and  simulated  ocean  salinity  and  tem- 
perature structure  for  the  year. 
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I.     INTRODUCTION 

Salinity  in  the  oceanic  mixed  layer  and  its  change  in  depth  with  time  can  have  a 
significant  impact  on  upper  ocean  dynamics.  Alterations  of  upper  ocean  salinity  caused 
by  precipitation  and  evaporation  through  buoyancy  fluxes  may  also  change  the  tem- 
perature profile  of  the  ocean,  altering  the  sound  velocity.  The  purpose  of  this  study  was 
to  examine  the  effects  of  discrete  precipitation  events  on  the  evolution  of  the  mixed  layer 
depth  (MLD)  at  Ocean  Station  "P"  (50°N,  145°W). 

Previous  one-dimensional  mixed  layer  prediction  experiments  had  failed  to  accu- 
rately explain  and  predict  the  behavior  of  salinity  with  time  and  depth.  Possible  reasons 
for  this  were: 

1.  the  amount  of  precipitation  and  evaporation  was  incorrectly  specified, 

2.  the  precipitation  and  evaporation  distribution  with  time,  relative  to  wind  events 
was  incorrectly  specified,  and  or 

3.  the  advection  was  incorrectly  specified  or  neglected. 

This  study  focused  specifically  on  the  variation  of  precipitation  and  evaporation  with 
time.  The  total  amounts  of  precipitation  and  evaporation  are  also  varied.  Possible  ef- 
fects of  advection  are  discussed. 

A.     LITERATURE  REVIEW 

Kraus  and  Turner  (1967)  provided  the  fundamental  turbulent  kinetic  energy  (TKE) 
theory  on  which  many  subsequent  oceanic  mixed  layer  models  are  based.  This  model 
considered  the  turbulent  kinetic  energy  equation  in  a  one-dimensional  mixed  layer 
scheme,  using  the  approximately  decoupled  state  of  the  mechanical  and  thermal  energy 
equations.  It  was  assumed  that  the  heat  input  at  the  air-sea  interface  and  the  mass 
entrainment  at  the  base  of  the  mixed  layer  were  uniform  throughout  the  mixed  layer  for 
short  time  scales.  Wind  stress  effects  on  the  air-sea  interface  acted  to  increase  the  po- 
tential energy  in  the  mixed  layer  by  generating  turbulent  kinetic  energy  for  mixing  and 
entrainment.  Entrainment  was  considered  to  be  proportional  to  the  wind  energy  input 
to  the  water  column  minus  the  work  performed  in  overcoming  the  buoyant  forces 
throughout  the  mixed  layer.  This  type  of  turbulent  erosion  modeling  considered 
buoyancy  and  wind  forcing  only,  and  it  neglected  viscous    dissipation.  This  neglect  of 


dissipation  has  generally  been  considered  a  serious  flaw.  Kraus  and  Turner  (1967)  did 
not  consider  salinity  effects  in  their  mixed  layer  model. 

Subsequent  works,  such  as  Geisler  and  Kraus  (1969),  Miropol'skiy  (1970),  and 
Denman  (1973)  have  attempted  to  consider  dissipation  by  assuming  dissipation  to  be  a 
fixed  portion  of  wind  stress.  The  Denman  (1973)  model  was  used  in  the  first  mixed  layer 
simulation  experiment  to  use  Ocean  Station  "P"  data  (Denman  and  Miyake,  1973). 

Miller  (1976)  was  the  first  to  include  salinity  and  its  effects  on  the  density  structure 
in  a  mixed  layer  model.  Miller  modified  Denman  (1973)  to  include  salinity  and  used  data 
from  BOMEX  Period  III  upper  ocean  soundings  for  periods  of  as  long  as  10  days  for 
forcing  and  initial  conditions.  Four  specific  cases  were  examined  by  selecting  distinctly 
different  initial  conditions  for  the  salinity  and  temperature  profiles  and  forcing  data  from 
BOMEX  surface  observations.  In  three  cases,  Miller  examined  the  effects  of  three  dif- 
ferent initial  conditions  while  the  fourth  case  included  a  precipitation  event. 

The  Miller  study  indicated  that  mixed  layer  salinity  greater  than  the  salinity  of  the 
water  below  the  mixed  layer  may  cause  convective  overturning  due  to  density  instability 
unless  a  sufficiently  large  temperature  decrease  occurred  simultaneously  across  the 
interface.  A  large  density  jump  across  the  interface  acted  to  slow  the  mixed  layer 
deepening  rate  due  to  an  increase  in  potential  energy.  Without  entrainment  salinity  flux 
effects,  the  cooling  induced  by  entrainment  was  greater  than  the  heating  from  the  sur- 
face, resulting  in  net  mixed  layer  cooling.  With  salinity  included  in  the  model,  the  layer 
deepened  more  slowly  and  entrained  less  cold  water  at  the  bottom,  allowing  surface 
heating  to  dominate  and,  therefore,  the  mixed  layer  temperature  to  increase.  Overall,  the 
results  of  these  cases  demonstrated  the  importance  of  salinity  fluxes  and  the  need  to  in- 
clude the  salinity  profile  in  the  prediction  of  mixed  layer  deepening  and  resultant  tem- 
perature structure. 

Miller  demonstrated  the  potential  significance  of  precipitation  episodes  in  mixed 
layer  modeling  in  a  fourth  case.  His  study  demonstrated  that  precipitation  can  contrib- 
ute to  salinity  fluxes  within  the  mixed  layer.  However,  the  fact  that  salinity  regularly 
fluctuates  by  ±0.15%  (due  to  horizontal  advection)  makes  decreases  in  salinity  due  to 
rainfall  difficult  to  detect  in  observations.  Though  Miller  admittedly  did  not  understand 
the  nature  of  the  effects  of  rainfall  on  the  mixed  layer,  he  did  propose  that  even  minor 
changes  in  the  salinity  flux  due  to  precipitation  may  have  a  long  term  effect  on  the  mixed 
laver  thermal  structure. 


Paulus  (1978)  examined  the  effects  of  a  salinity  profile  on  density  structure  using  a 
modified  version  of  the  Camp  (1976)  one-dimensional  mixed  layer  model.  The  effects 
of  a  surface  salinity  flux  on  the  mixed  layer  were  examined  using  a  constant  evaporation 
and  precipitation  rate.  Model  initialization  was  accomplished  with  hydrocast  data  or 
XBT  data  with  historical  salinity  profiles  adjusted  to  realistically  correspond  to  the 
temperature  profiles.  The  incorporation  of  observed  salinity  profiles  into  the  model  had 
little  effect  on  sea  surface  temperature  (SST)  and  MLD  during  the  summer  regime  (when 
MLD  was  shallow).  However,  when  using  the  same  forcing  conditions  for  the  winter 
regime  the  inclusion  of  a  salinity  profile  acted  to  produce  a  warmer,  shallower  mixed 
layer  than  obtained  when  density  was  only  a  function  of  temperature. 

Paulus  suggested  that  the  retreat  of  the  isothermal  layer  could  be  enhanced  due  to 
precipitation  by  slowing  the  deepening  rate  and  shoaling  tendency  of  the  mixed  layer 
under  downward  heat  flux.  With  respect  to  salinity  flux,  Paulus  concluded  that  imbal- 
ances in  evaporation  and  precipitation  had  more  pronounced  effects  on  SST  and  MLD 
in  the  long  term  than  in  the  short  term.  His  study  concluded  that  this  was  a  cumulative 
effect.  Therefore,  Paulus  concluded  that  the  addition  of  salinity  effects  to  the  mixed 
layer  model  was  desirable. 

Garwood  (1977)  introduced  a  one-dimensional  bulk  model  of  the  oceanic  mixed 
layer  that  included  dynamic  instability  with  turbulent  erosion.  This  model  modulated 
the  mean  entrainment  rate  by  diurnal  heating  and  limited  maximum  layer  depth  by  en- 
hancing the  effect  of  viscous  dissipation.  This  allowed  a  cyclic  steady  state  to  occur  over 
an  annual  time  period.  It  recognized  the  important  effects  of  long  term  salinity  fluxes 
on  the  mixed  layer  thermal  profile  due  to  surface  heat  fluxes,  precipitation,  and 
entrainment  at  the  base  of  the  mixed  layer.  This  was  the  model  of  choice  for  this  study. 

Garzon  (1987)  used  an  oceanic  mixed  layer  model  to  determine  the  specific  effects 
of  fresh  water  flux  on  the  seasonal  pycnocline  in  the  mixed  layer  at  Ocean  Station  "P." 
First,  sensitivity  of  the  model  proposed  by  Garwood  (1977)  was  examined  by  performing 
a  hypothetical  experiment  using  constant  forcing  with  a  wide  range  of  precipitation  and 
evaporation  values.  These  results  were  compared  to  an  identical  model  run  with  P-E 
(precipitation  minus  evaporation)  =  0.0.  This  simple  hypothetical  case  examined  the 
relative  importance  of  excess  evaporation  or  precipitation  over  a  time  period  of  one- 
year.  For  the  case  with  excess  precipitation,  a  reduction  in  downward  buoyancy  and  a 
decrease  in  surface  salinity  occurred,  resulting  in  a  shallower  MLD  than  would  be  ex- 
pected if  P-E  =  0.     For  the  case  with  excess  evaporation,  an  increase  in  downward 


buoyancy  flux  and  an  increase  in  surface  salinity  occurred.    This  yielded  a  deeper  MLD 
than  would  occur  with  P-E=  0. 

For  the  annual  simulation  at  Ocean  Station  "P,"  Garzon,  using  data  from  1967, 
modeled  temperature  and  salinity  profiles  with  P-E  =  0.0m  year,  P-E  =  0.5m  year,  and  for 
monthly  representative  values  of  precipitation.  Overall,  MLD  and  surface  salinity  from 
the  model  well  approximated  those  of  the  actual  data.  However,  near  the  end  of  the 
study  period,  MLD  and  salinity  values  demonstrated  increasing  error  particularly  in  ap- 
proximating significant  anomalous  "spikes."  Garzon  suggested  that  a  more  realistic 
representation  of  discrete  precipitation  and  evaporation  events  would  result  in  more  ac- 
curate model  output,  especially  for  the  long  term.  Additionally,  Garzon  adjusted  model 
output  for  advection  of  salinity  by  subtracting  the  linear  fit  from  the  predicted  salinity 
time  series. 

Precipitation  and  evaporation  events  have  been  demonstrated  to  have  a  significant 
impact  on  mixed  layer  depth  on  both  short  and  long  time  scales.  These  factors  can  act 
alone  or  may  be  coupled  with  such  events  as  synoptic  scale  storms  to  modify  the  effects 
of  these  events.  On  a  short  time  scale,  daytime  heating  tended  to  cause  shallowing  of 
the  oceanic  mixed  layer  due  to  heating  of  the  mixed  layer  despite  the  occurrence  of 
evaporation  which  tended  to  promote  deepening.  Upward  heat  flux  broke  down  the 
mixed  layer  causing  deepening.  Precipitation  during  this  time  period  promoted  shal- 
lowing, thus  opposing  this  effect.  These  diurnal  changes  in  mixed  layer  depth  were 
moderated  by  annual  effects  that  were  associated  with  the  annual  radiation  cycle. 
Strong  winds,  often  associated  with  large  scale  synoptic  events,  when  combined  with  low 
insolation  values  in  winter,  yielded  deep  mixed  layer  depths.  However,  precipitation 
events  often  accompanying  these  winds  acted  to  lessen  the  effects  of  winds  and  were 
associated  with  discrete  shallowing  episodes  that  last  for  short  periods  of  time.  These 
effects  were  cumulative.  Just  as  a  reduction  or  increase  in  wind  speed  in  winter  yielded 
a  substantially  shallower  or  deeper  mixed  layer  depth  in  summer,  a  reduction  or  increase 
of  rain  during  the  same  period  yielded  a  deeper  or  shallower  mixed  layer  depth  in  sum- 
mer. 

B.     RESEARCH  PLAN:  LOCATION  AND  PERIOD  OF  STUDY 

For  this  study,  a  time  period  was  chosen  to  provide  the  best  overlap  of  available 
data  bases.  The  starting  date  of  23  November  196S  with  a  study  period  of  400-days  was 
selected  for  this  study,  since  it  provided  data  base  overlap,  that  is,  surface  forcing  and 


CTD  data  was  available  for  this  time  period.  This  numerical  simulation  of  400-days  used 
a  time  step  of  one  hour. 

The  study  area  was  between  49°  N  and  51°  X,  and  144°  W  and  146°  W  centered  at 
Ocean  Station  "P"  (Fig.  1  on  page  6).  From  1949  to  1981  the  waters  around  Ocean 
Station  "P"  have  been  occupied  by  ocean  station  vessels  year-round.  This  provided  nu- 
merous data  bases,  such  as,  XBT,  CTD,  and  atmospheric  observations  for  analysis. 

Tabata  (1961,1965)  described  oceanographic  and  meteorological  conditions  at  Sta- 
tion "P"  in  detail.  The  dominate  pressure  systems  located  in  the  North  Pacific  Ocean 
are  the  low  pressure  system  or  Aleutian  Low  in  winter  and  the  high  pressure  system  or 
North  Pacific  High  in  summer.  These  two  systems  govern  the  distribution  of  mean 
surface  winds  over  the  Northeast  Pacific  Ocean.  In  winter  the  Aleutian  Low  directs  air 
flow  northward  into  the  Gulf  of  Alaska,  while  in  summer  the  North  Pacific  High  directs 
winds  southeast  at  Ocean  Station  "P."  Additionally,  monthly  mean  wind  speeds  are 
approximately  twice  as  large  in  magnitude  (13ms)  in  winter  as  in  summer  (6m  s). 

Station  "P"  is  located  in  a  distinct  oceanic  regime  (Fig.  2  on  page  7).  The  surface 
waters  of  the  region  are  under  the  influence  of  the  weak  westward  flow.  Specifically,  this 
flow  is  composed  of  the  West  Wind  Drift  and  the  recirculated  waters  from  the  Alaska 
Stream.  The  former  carries  mixed  waters  of  the  Kuroshio  and  the  Oyashio,  which  are 
generally  indistinguishable  by  the  time  that  they  reach  Ocean  Station  "P."  To  the  east 
of  this  area,  the  current  turns  north,  forming  the  Alaskan  current. 

The  vertical  structure  of  this  water  is  characterized  by  three  zones:  (1)  the  upper 
"seasonal"  zone  (0-100m),  (2)  the  halocline  (100-200m),  and  (3)  the  lower  zone  (>  200m). 
Within  the  upper  zone  the  salinity  is  relatively  low  (~  32.7ppm).  With  respect  to  sea- 
sonal effects,  in  winter  the  temperature  and  salinity  in  zone  1  are  relatively  homogeneous 
and  isothermal.  In  summer  a  significant  vertical  salinity  gradient  exists  in  this  upper- 
most zone,  as  well  as  a  well  defined  thermocline,  with  a  temperature  decrease  with  depth 
of  as  much  as  8°C  per  20m.  In  the  halocline,  zone  2,  salinity  increases  with  depth  by 
as  much  as  l.Oppm  per  100m.  In  zone  2  temperature  decreases  with  depth  typically, 
unless  an  inversion  is  present.  In  the  lower  zone  salinity  increases  only  gradually  with 
depth,  while  temperature  decreases  gradually  with  depth  (Fig.  3  on  page  8). 

Tabata  (1961)  indicates  that  there  is  an  annual  average  excess  of  precipitation  over 
evaporation  of  ~  0.5m.  Garzon  (1987)  also  estimates  an  excess  annual  precipitation. 
Despite  this  annual  excess  of  precipitation,  salinity  within  the  water  column  does  not 
decrease    on    an  annual  time  scale,  although    it    can  vary  as    noted  by  Tabata  (1961). 
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Fig.   1.      Region  Surrounding  Ocean  Station  "P"  (from  Tabata,  1961) 


Fig.  2.      Oceanic  Regime  at  Station  "P."  (from  Favorite,  1976) 


Fig.  3.      Example  Temperature  and  Salinity  Profiles 


especially  within  the  more  dynamic  mixed  layer.  Horizontal  advection  must  account  for 
this  readjustment.  As  noted  by  Tabata,  upwelling  is  weak  and  nearly  insignificant  within 
these  waters,  when  compared  with  horizontal  advection.  As  a  result  this  salinity  balance 
appears  to  be  maintained  by  horizontal  advection  as  more  saline  water  is  transported 
into  the  area  by  ocean  currents. 

Despite  the  consistent  westerly  current  coming  into  this  area,  anomalous  events 
have  occurred  resulting  in  horizontal  transport  of  water  into  the  area  from  a  non- 
westerly  direction.  Tabata  (1961)  noted  an  unusual  increase  in  temperature  and  decrease 
in  salinity  within  the  upper  zone  and  halocline  during  spring  1957  through  summer  1958. 
This  was  attributed  to  transport  in  water  lying  south  of  Ocean  Station  "P"  into  the  study- 
area.  The  suggested  cause  for  this  transport  was  the  shift  of  winds  in  the  region  associ- 
ated with  unusual  location  of  the  atmospheric  pressure  system  over  the  northeast  Pacific 
Ocean.  This  displacement  of  the  North  Pacific  High  yielded  westerly  instead  of  the  usual 
southeasterly  winds,  which  could  have  generated  southerly  drift  currents  in  the  vicinity 
of  the  station. 

Overall,  the  currents  in  this  area  are  weak.  Except  for  occasional  anomalous  trans- 
port associated  with  unusual  synoptic  weather  systems,  salinity  follows  a  roughly  annual 
cycle  maintained  by  horizontal  advection.  and  the  waters  of  the  study  area  are  relatively 
homogeneous  horizontally.  These  factors  combine  to  provide  an  excellent  region  for  a 
one-dimensional  modeling  study  of  the  thermodynamics  of  the  upper  ocean. 


II.     DATA  SOURCES  AND  PREPARATION 

A.  DATA  SOURCES 

This  study  was  performed  using  various  sources.  Atmospheric  forcing  data  was  or- 
iginally obtained  from  the  National  Weather  Records  Center,  Asheville,  North  Carolina. 
This  data  included  three-hourly  observations  of  wind  direction  and  speed,  air  temper- 
ature, sea  surface  temperature,  dew  point  temperature  and  cloud  coverage  in  octaves. 
Atmospheric  data  was  interpolated  from  three  hourly  to  one  hourly  data  for  computa- 
tion of  hourly  surface  fluxes.  Precipitation  forcing  data  was  not  directly  available  and 
was  derived  as  explained  later  in  this  chapter. 

Two  types  of  oceanic  data  were  used.  CTD/STD  data  was  provided  by  the  Marine 
Environmental  Data  Service,  Ottawa,  Canada.  These  salinity  and  temperature  casts 
were  digitized  from  depths  of  0  to  1500  m  at  various  intervals,  with  more  concentrated 
data  provided  within  the  mixed  layer.  These  casts  were  not  taken  at  regular  intervals, 
but  were  taken  as  operationally  feasible.  There  were  119  casts  available  for  analysis 
during  the  study  period  from  station  "P."  Bathythermograph  casts  were  not  available 
for  this  time  period. 

Three  hourly  atmospheric  forcing  data  was  used  in  this  study  to  compute  the  fluxes 
of  heat,  radiation  and  momentum  necessary  for  forcing  the  mixed  layer  model  with  an 
hourly  time  step.  Specifically,  these  boundary  conditions  were  wind  stress  in  two- 
dimensional  components,  incident  solar  radiation,  heat  loss  (long  and  short  wave),  and 
precipitation  and  evaporation.  The  calculation  of  precipitation  will  be  discussed  later. 
Temperature  and  salinity  profiles  from  STD  taken  23  November  1969  were  used  to  ini- 
tialize the  model  profiles,  and  CTD  STD  casts  were  used  for  verification.  Bucket  tem- 
peratures were  used  for  comparison  of  with  predicted  temperature. 

B.  VERTICAL  INTERPOLATION 

The  oceanic  profiles,  from  both  CTD/STD  sources  were  vertically  interpolated  from 
0  to  199  m  to  a  standard  grid  of  1  m  spacing.  Mixed  layer  depth,  mixed  layer  temper- 
ature, and  mixed  layer  salinity  were  derived  from  these  interpolated  profiles.  There  is 
not  a  standard  way  to  compute  mixed  layer  depth  from  a  CTD/STD  profile.  For  this 
study  the  mixed  layer  depth  has  been  defined  to  be  that  depth  at  which  the  temperature 
profile  changes  exceed  0.10°  C  in  5  m.   Mixed  layer  temperature  and  mixed  layer  salinity 
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values  were  obtained  by  averaging  temperature  from  1  m  to  the  mixed  layer  depth  as 
determined  above. 
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III.     THEORY 

A.  OCEANIC  MIXED  LAYER  STRUCTURE 

A  familiarization  with  the  definitions  associated  with  the  oceanic  mixed  layer  ther- 
mal structure  is  required  to  understand  the  dynamics  involved  in  mixed  layer  modeling. 
As  defined  in  Garzon  (1987)  the  oceanic  mixed  boundary  layer  is  assumed  to  be  the  fully 
turbulent  region  of  the  upper  ocean  bounded  above  by  the  air-sea  interface  and  below 
by  a  relatively  static  water  mass.  The  mixed  layer  is  considered  to  be  "quasi- 
homogeneous,"  that  is,  it  is  approximated  to  be  well  mixed.  Thus  temperature  within 
the  mixed  layer  is  considered  to  be  relatively  constant.  Fig.  4  on  page  13  graphically 
depicts  an  oceanic  vertical  temperature  profile.  The  turbulent  mixed  layer  of  depth  h  is 
contained  within  the  upper  most  oceanic  region,  below  which  lies  the  turbulent 
entrainment  zone,  from  depths  z=-h  to  z=-h-d>,  with  z  positive  upward.  Within  the 
entrainment  zone  there  is  a  temperature  jump  of  AT .  Below  the  entrainment  zone  is  the 
quasi-stable  water  mass  and  here  temperature  decreases  with  increasing  depth. 

Mixed  layer  structure  (Fig.  5  on  page  14)  for  salinity  is  similar  to  that  for  temper- 
ature (Fig.  4  on  page  13).  There  is  a  relatively  constant  salinity  within  the  turbulent 
mixed  layer  from  z=0  to  z=-h,  a  salinity  jump  AS  across  the  entrainment  zone  from 
z=-h  to  z  =  -h-S,  and  a  gradually  increasing  salinity  with  depth  below  the  mixed  layer. 

An  oceanic  mixed  layer  prediction  scheme  requires  an  accurate  (or  near  accurate) 
assessment  of  the  deepening  or  shallowing  rate  of  that  layer.  Garwood  (1977)  accounted 
for  both  of  these  processes.  The  equations  used  by  Garwood  to  best  describe  these 
processes  are  developed  as  follows. 

B.  ENTRAINMENT 

The  Garwood  (1977)  one-dimensional  mixed  layer  model  considers  both  turbulent 
erosion  and  dynamic  instability.  Fundamental  to  the  solution  of  this  model  is  an  accu- 
rate assessment  of  the  entrainment  velocity,  IV,.   Starting  with: 

^L=We-W{-h)  (3.1) 

Assuming  no  mean  vertical  motion,  W  (-h)  =  0,  this  equation  becomes: 

»>■ fr  (3-2) 

CI 
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Fig.  4.      Example  Temperature  Profile 

The  solution  for  W,  can  be  found  by  solving  the  total  turbulent  kinetic  energy  (TKE) 
budget  written  as: 


d    E 
dt    2 

_v,w,c±_ 

cz 

cz 

H'griv]-^ 

_,f"(-f^r)]-£(3-3) 

where: 
E  =  JT'2  +  V'i  +  W\  is  the  Total  TKE. 

(see  Table  1  on  page  20  for  definitions  of  constants) 

The  five  terms  in  Equation  (3.3)  from  left  to  right  are: 
1.  the  time  rate  of  change  of  turbulent  kinetic  energy, 
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Fig.  5.      Example  Salinity  Profile 

2.  shear  production, 

3.  buoyancy  flux, 

4.  turbulent  diffusion  (or  redistribution),  and 

5.  viscous  dissipation. 

This  equation  may  be  examined  for  the  entrainment  zone.  Within  the  entrainment  zone, 
dissipation  and  shear  production  are  assumed  negligible.  This  produces  a  mean  turbu- 
lent kinetic  energy  equation  for  the  entrainment  zone  of  the  form 


JLJL  =  agTiy-JLiy 

ct    2  cz 


2       Po 


(3.4) 
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across  the  entrainment  zone  of  thickness  <5.    Integrating  equation  (3.4)  yields: 


-h-6  J-/l"<5  J 


(3.5) 


Assuming  steady  state,  the  time  rate  of  change  of  turbulent  energy  may  be  neglected. 
Equation  (3.5)  becomes: 


Jw^^-Rt+tFJ] 


f*£YL- 


(3.6) 


Solving  this  integral  produces  equation  (3.7) 
agS 


{rii'\_h  +  riv'\_h_6)-iv'(f  +  -^y-_l6  =  o  (3.7) 

When  specifically  examining  this  equation  at  the  base  of  the  mixed  layer  (or  the  top  of 
the  entrainment  zone): 


c    „„,  E   ,    F 


,griv\_^f7ir{f  +  ^)\_ 


(3.8) 


Kraus  and  Turner  simplified  the  solution  at  this  point  by  neglecting  the  redistribution 
term  and  combining  the  entrainment  mean  kinetic  energy  equation  with  the  steady  state 
turbulent  kinetic  energy  equation  for  the  mixed  layer.  Garwood  does  not  neglect  this 
term,  but  estimates  its  magnitude  with  the  following  parameterization: 


J-iy.fJL+PL^  <E> 

cz  "  (  2  +  Po  )l-h~      re 


(3.9) 


Thus  the  redistribution  term  is  equal  to  the  mean  turbulent  energy  available  near  the 
base  of  the  mixed  layer  divided  by  the  time  scale  required  to  transport  this  energy  into 
the  area.   Parameterizine  this  time  scale  as 


(3.10) 


W 


where  |  W  I  is  the  root-mean  square  vertical  turbulent  velocity,  |  W  \  =  (W2)1'2,  and  the 
time  scale  is  approximated  as  the  absolute  value  of  the  mixed  layer  depth  divided  by  the 
absolute  value  of  the  vertical  velocitv. 
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Using  Equations  (3.8  to  3.10)  and  assuming  that  TW  |_*  =  ATWt , 

g(a.AT-  PAS)l)  e  = 1 


(3.11) 


Here  vertical  transport  of  turbulence  exactly  balances  entrainment  buoyancy  flux.   This 
equation  may  be  solved  for  the  entrainment  velocity: 


I  W  1  <  E  > 

e      hg(aAT-pAS) 


(3.12) 


To  solve  this  equation  for  W,  the  values  of  <  E  >,  \W'\  ,  T,  h,  and  AT  must  be 
known  while  g  and  a  are  constants  (Table  1  on  page  20).  From  Garwood  (1977)  the 
total  turbulent  kinetic  energy  is  solved  by  using: 


ct 


(h<E>)=GQ+Gh-B0-Bh-D 


(3.13) 


In  order  to  solve  equation  (3.13),  the  following  equations  for  shear  production  at 
z  =  0  and  z=-h,  and  buoyancy  damping  production  at  z  =  0  and  z=-h,  and  dissipation 
are: 

wind  shear  production, 


G0  =12 


$y 


■^(f 


3/4 


entrainment  shear  production, 


Gh  =  [  <  u  >2  +  <  v  >2]  We, 
where  the  x  and  y  components  of  the  mixed  layer  current  are: 


1.  —  {h  <  u  >  )  =fh  <  v  >  +  -r-,  and 

dt  P 

2.  ■j^(h<v>)=Jh<u>  +  ^r. 


(3-14) 


(3.15) 
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buoyant  damping/production  due  to  surface  buoyancy, 


Bn  =  a° 


<2c 


pcP 


-PgS(E-P), 


(3.16) 


buoyant  damping  due  to  entrainment, 


Bh  =  aghATlVe,  and 


(3.17) 


dissipation, 


D=2<E> 


3/2 


(3.18) 


C.     MIXED  LAYER  THERMODYNAMICS:  TEMPERATURE  AND  SALINITY 
BUDGETS 

Wind  stress  generated  turbulence  at  the  air-sea  interface  provides  the  energy  to  de- 
stabilize and  entrain  the  water  below.  This  mixing  process  deepens  the  mixed  layer  and 
mixes  temperature  and  salinity  between  the  surface  and  the  entrainment  zone.  Thus,  the 
conservation  of  heat  requires  that: 


cT 
8t 


cz 


(3.19) 


Vertically  integrating   Equation  (3.19)  across  the  entrainment  zone,  z  =  —  h  —  S   to 
z  =  —  h,  and  assuming  T  W  vanishes  below  z  =  —  h  —  5  gives 


-ATWe=T'lV'{-h) 


(3.20) 


where: 
AT  =  T  (z  =  -  h)  -  T(z  =  -h-  5). 


Vertically  integrating  Equation  (3.19)  across  the  entire  mixed  layer,  z=-h  to  z  =  0 


gives 


,o     _ 


cT    ,        r          8TW    j 
—r-dz=  - — dz 

-h  0[  J-h~6         oz 


(3.21) 


Using  Equation  (3.20)  and  assuming  T  to  be  well  mixed  in  z,  Equation  (3.21)  becomes 
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dT     riv(Q)      riv(-h) 


ct  h  h 


+ i (3.22) 


or 


ct        pCph 

ATWe 
h 

Here  : 

1.  &T=T- 

-T{-h-  S),  and 

(3.23) 


2-  Q0lpCp  =  -TlV'(0). 
The  one-dimensional  salinity  budget  for  the  oceanic  mixed  layer  is: 


cS  61VS' 


ct  cz 


(3.24) 


It  is  assumed  that: 

1.  S(z)  is  well  mixed  from  the  bottom  of  the  mixed  layer  to  the  top  of  the  mixed  layer, 

2.  S(z)  has  a  near  discontinuity  between  z=  -h  and  z=  -h-<5,  and 


3.    the  vertical  flux  of  S  (or  S'  W)  vanishes  at  z  <  —h  —  b. 

Fig.  5  on  page  14  depicts  this  condition.  Here  the  turbulent  region  within  the  mixed 
layer  yields  a  near  homogeneous  state  with  respect  to  salinity.  Examining  Equation 
(3.24)  and  integrating  gives: 


4~dz  = 


-h(i)-6   CI 


0  \        dz  (3.25) 

-h-6       0Z 


With  h  variable  in  time,  application  of  the  Leibnitz  Rule  of  integration  and  assuming  <5 
is  vanishing  small, 

tt  dh 


AS  -TT-  =  -  IV'S'(  -h)  +  IV'S'(  -h-5).  (3.26) 


Assuming  that  turbulence  is  negligible  below  the  mixed  layer,  i.e.  W'S'{  —  h  —  S)  ~  0. 
This  makes  the  "jump  condition"  for  salinity  entrainment  into  the  mixed  layer: 

~  dh 


AS-±r-  =  -lV'S'{-h).  (3.27) 
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where: 
AS  =  S(z  =  -  h)  -  S{z  =  -  h  -  S) 

D.      SALINITY  EFFECT  ON  MIXED  LAYER  DYNAMICS 

1.     The  Shallowing  Effect 

Kraus  and  Turner's  (1967)  model  provides  an  approximation  for  the  non- 
entraining  or  shallowing  mixed  layer,  We  =  0, 

h=1£l b. =L  (3.28) 

Q    ugQol{pCp) 

Here  L  is  the  Obokov  length  scale  and  equation  (3.28)  is  the  solution  when  there  is  no 
entrainment.  The  surface  buoyancy  flux  is  B0  =  {agq0l(pCp))  and  U.  is  the  surface  friction 
velocity,  with  LI  =  (pJp)CDU2  .  The  value  of  h  is  inversely  proportional  to  Q0.  The 
mixed  layer  temperature  is  even  more  sensitive  to  Q0  than  is  h: 

4f  =  Qol(pCph)  ozQ2oc-\  (3.29) 

If  the  buoyancy  flux  attributable  to  precipitation  and  evaporation  is  included,  Equation 
(3.28)  becomes 


C  ul 


Q    agQ0l{pCp)  +  fig{P  -  E)S 


(3.30) 


Where  the  salinity  concentration  coefficient,  /?  =  0.76x1 0-3Kglgm  ,  and  P-E  is  the  net 
precipitation  minus  evaporation  in  msec.  With  the  salinity  flux  included,  the  surface 
buoyancy  flux  becomes: 

#o  =  «gQoiPCp  +  Pg(P  -  E)S,  (3.31) 

and  the  effect  on  mixed  layer  temperature  becomes: 


II 

dt 


HT  -  (QolpCp)lhaQ0(Q0  +  fila{P  -  E)S).  (3.32) 


As  a  result,  precipitation  minus  evaporation  influences  mixed  layer  temperature  indi- 
rectly. For  positive  Q0  (net  heat  into  the  surface),  a  positive  (P-E)  will  decrease  h,  con- 
centrating the  heat  into  a  shallow  layer.  For  the  case  where  evaporation  exceeds 
precipitation  a  decrease  of  mixed  layer  temperature  with  respect  to  time  results. 
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2.     The  Deepening  Effect 

For  entraining,  deepening  mixed  layers,  the  buoyancy  discontinuity  at  the  bot- 
tom of  the  mixed  layer  must  be  augmented  to  include  a  possible  salinity  jump: 


where: 


Ab  =  ag&T -  pgAS 


AS  =  S-S(-h-6), 

S  is  the  mixed  layer  salinity,  and 

S(  —  h  —  S)  is  the  salinity  just  below  the  mixed  layer. 
The  mixed  laver  salinitv  budset  is: 


(3.33) 


cS 

01 


S(E-P) 

h 


WeAS 


(3.34) 


Equation  (3.34)  demonstrates  that  the  change  in  salinity  with  respect  to  time 
depends  on  fresh  water  flux  at  the  surface,  entrainment  velocity,  salinity  advection, 
mixed  layer  depth,  and  upon  the  salinity  jump  at  the  base  of  the  mixed  layer. 

Table   1.      MODEL  CONSTANTS 


Model  Constants 

constant 

definition 

value 

a 

thermal  expansion  coefficient 

-  W  :Kl 

/? 

salinity  expansion  coefficient 

~  8  x  lO^ppnr1 

c, 

tuning  coefficient 

-  1 

c 

tuning  coefficient 

^  .2 

cF 

specific  heat 

~  40007  •  kg-'  .  A'"1 

p 

sea  water  density 

-  1026.5%.  m~3 
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IV.     HYPOTHETICAL  EXPERIMENTS 

A.     EFFECTS  OF  DISCRETE  PRECIPITATION  EVENTS 

The  purpose  of  this  study  was  to  determine  the  lasting  effects  of  discrete  precipi- 
tation events  on  the  mixed  layer.  For  "case  1"  in  this  study  short-term  precipitation 
events  of  10  cm  in  one  day  with  no  evaporation  was  applied  monthly  to  the  Garwood 
(1977)  model  over  a  365  day  study  period.  This  yielded  twelve  separate  events  for  de- 
tailed study.  At  all  other  times  during  the  year  both  evaporation  and  precipitation  were 
assumed  to  zero.  The  resulting  output  from  the  model  was  examined  individually  for 
30  day  periods  following  each  rain  event  to  better  reveal  long  term  trends.  The  precipi- 
tation event  was  applied  on  the  fourth  day  of  each  monthly  period.  Actual  temperature 
and  salinity  profiles  taken  from  a  CTD  cast  for  November  23,  1968  were  used  as  initial 
conditions.  Heat  and  momentum  fluxes  were  computed  from  the  actual  surface 
meteorological  and  observations. 

Fig.  6  on  page  22  and  Fig.  7  on  page  23  are  a  conglomeration  of  the  results  of  these 
twelve  separate  events.  Other  figures  in  this  chapter  represent  the  individual  events  that 
were  selected  for  closer  examination  at  a  60  day  period. 

The  effects  of  precipitation  events  on  the  MLD  are  graphically  represented  by  the 
figures  in  this  chapter.  The  equations  developed  in  Chapter  3  will  be  used  to  explain  the 
results  and  are  repeated  here: 


6T         Qo         AFW 


/ 


e 


ct        pCJi 


(4.1) 


The  change  in  mixed  layer  temperature  (T)  with  respect  to  time  is  a  function  of: 

1.  MLD,  h, 

2.  entrainment  rate,  Wt 

3.  temperature  jump  at  the  base  of  the  mixed  layer,  AT,  and 

4.  net  downward  heat  flux,  Q0. 

For  a  case  of  no  entrainment,  Wt  =  0,  the  temperature  will  rise  for  Q0  >  0,  and  it  will  fall 
for  <2o  <  0.  The  rate  of  rise  and  fall  is  inversely  proportional  to  h.  If  Q0  <  0  and  there 
is  entrainment,  Wt  >  0,  the  mixed  layer  temperature  can  only  decrease  if  AT>0.    The 


21 


Eo 
o 
—  ^^  r\] 


«3  ° 

1   CM 


v        r 


b. 


&5  o 


i        I        I 1 1 1 1 — 

25     50     75     100    125    150    175 


£ 


-K- 


h7  o 


C'  §,° 


25 


50 


75 


100 


in 


i 


Y 


Go   o 


25 


50 


75 


100 


a.  -5-  ° 


',^^/ifiir'ffffiWii^,,,, 


0' 


25 


50 


75 


At 


fifUj 


100 


i     r r 

125    150    175 


e. 


CO  C\] 


ID 


o 


'J 


25     50     75     100    125    150    175 


f.    ° 


a 


o 


8 

o 
o 


r^^AiVrAV^ 


200 


125    150    175    200 


125    150    175    200 


200 


75     100    125    150    175    200 

TIME(DY) 


~      I      I      I 1 1 1 1 

25     50     75     100    125    150    175    200 


Fig.  6.  Conglomeration  of  SL\  Events  (Dec  1968  to  May  1969):  (a)  differential 
MLD,  (b)  dillercntial  temperature,  (c)  differential  salinity,  (d)  MLD  with 
L>P  =  U,  (e)  wind  speed,  and  (I)  net  surface  heat  flux. 
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Fig.  7.      Conglomeration  of  Six  Events  (Juii  to  Nov  1969):     (a)  differential  MLD, 

(b)  differential  temperature,  (c)  differential  salinity,  (d)  MLD  with  E-P  =  0, 

(c)  wind  speed,  and  (0  net  surface  heat  (lux. 
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rate  of  decrease  is  greatest  for  small  h.   Various  factors  effect  \Ve,  and  will  be  considered 
later. 

Two  model  predictions  may  be  compared.  In  this  experiment  the  first  case  included 
the  discrete  precipitation  event,  while  the  second  case  had  no  events,  with  E-P  =  0  for  the 
entire  model  integration.  When  subtracting  the  results  of  one  study  using  equation  (4.1) 
from  another  such  study  for  the  same  time  using  the  same  equation,  the  results  provide: 

1.  differential  MLD  (A,  - /i2), 

2.  differential  temperature  (7",  -  T2),  and 

3.  differential  salinity  (Sx  —  S2). 


Differentiating  the  differential  temperature  with  respect  to  time  yields: 


c(T,  -  T2) 


ot 


h 


<2c 


Pc 


P        L 


Ar,H'el        AT2lVe2 


(4.2) 


Therefore,  in  order  to  provide  a  value  other  than  zero  for  8(T{  —  T2)jct  ,  the  differential 
temperature,  one  or  more  of  the  following  conditions  must  exist: 

1.  h^h2, 

2.  ireliH;,and 

3.  AT{^AT2. 

Again,  considering  the  many  variables  that  may  alter  the  entrainment  rate,  there  are 
numerous  factors  which  may  effect  this  differential  temperature,  making  interpretation 
of  even  a  simple  precipitation  event  complex,  but  more  easily  examined  than  the  case 
with  a  more  realistic  and  complex  E-P  time  series  of  forcing. 

Initially,  in  all  twelve  events  with  all  forcing  values  equal  (for  the  first  five  days),  the 
differential  temperature  ST=  Tx  —  T2  must  equal  zero.  This  can  be  seen  in  Fig.  6  on 
page  22  and  Fig.  7  on  page  23  for  the  first  five  days  of  each  30-day  period.  Upon  initi- 
ation of  the  precipitation  event  at  day  five  of  each  events  h{  =  h2  and  AT,  =  A7~2  ,  but 
W,  is  instantaneously  altered  in  the  precipitation  case.  In  this  case  with  E-P  <  0, there 
results  in  an  immediate  increase  of  surface  buoyancy,  B0  , 


B0  =  ag 


PCP 


-PgS(E-P). 


An  increase  in  B0  acts  to  decrease  the  total  TKE  from  equation  (3.13), 
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-|-(/z<£>)  =  G0  +  Gh-B0-  Bh-D 


and,  therefore,  entrainment  rate  from  equation  (3.12), 

W'\  <E> 


Wa  = 


gh{oAT -  pAS) 


As  We  decreases,  so  will  MLD,  which  acts  dynamically  to  shallow  the  MLD  relative 
to  the  case  with  P-E  =  0.  Thus  a  precipitation  event  increases  the  downward  surface 
buoyancy  through  decreasing  salinity,  which  ultimately  results  in  a  shallower  MLD. 
This  initial  result  should  be  apparent  throughout  the  twelve  experiments  where  there  is 
precipitation. 

While  the  immediate  change  in  MLD  through  reduced  entrainment  is  expected  on 
a  short  time  scale  (during  the  rain  event),  longer  term  changes  in  mixed  layer  dynamics 
may  result  from  a  single  short-term  precipitation  event.  In  response  to  mixed  layer 
shallowing,  MLD  and  the  temperature  jump  at  the  base  of  the  mixed  layer  will  change 
from  that  of  the  non-precipitation  case.  At  this  point  all  three  possible  conditions  re- 
quired to  yield  an  increase  in  the  differential  temperature  may  exist. 

Thus  analysis  of  the  effects  of  a  discrete  precipitation  event  requires  a  case  by  case 
examination.  Any  differences  in  MLD  will  yield  a  non-zero  magnitude  of  \hl  —  h2\ 
which  will  effect  the  differential  temperature.  A  fluctuation  in  differential  MLD  will  also 
moderate  the  effects  of  O0  ,  ATt  —  AT2,  and  IV el  —  Wt2  .  Any  difference  in  the  temper- 
ature jumps,  \ATX  —  AT2\  =£0  ,  will  also  yield  an  associated  change  in  the  differential 
temperature  and  a  transient  response,  as  well  as  the  altered  value  of  We . 

B.     ANALYSIS  OF  PRECIPITATION  EVENT  EFFECTS 
1.     Event  #  I  (Dec  1968) 

Fig.  6  on  page  22  represents  the  MLD  output  for  the  model  with  E-P  =  0.  This 
figure  indicates  that  this  event  began  in  the  late-fall  season  when  the  mixed  layer  was 
at  80  m  and  deepening.  Fig.  6  on  page  22  and  Fig.  8  on  page  26  indicates  that  there 
was  a  net  heat  loss  (Q0  <  0)  during  this  period  associated  with  a  mid-latitude,  winter  re- 
gime. There  were  moderate  winds  during  this  time  period  (U  =  5.0m/ s)  ,  and  maximum 
values  of  Um„  =  20.0m/ s  which  occurred  immediately  before,  during  and  after  initiation 
of  the  precipitation  event  of  day  5. 
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Fig.  8.      Event  #1  (Dec  1968):     (a)  difTerential  MLD  magnified  from  day  2  to  day 
8,  (b)  differential  temperature,  (c)  differential  salinity,  and  (d)  wind  speed 
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Prior  to  initiation  of  the  precipitation  event  all  differential  values  equaled  zero 
(Fig.  6  on  page  22  and  Fig.  7  on  page  23).  This  was  true  for  all  twelve  events.  Upon 
onset  of  precipitation  at  day  four,  MLD  decreased  from  65m  to  10m,  yielding  a  decrease 
in  the  differential  MLD.  This  was  coincident  with  shallowing  of  the  mixed  layer  asso- 
ciated with  a  lull  in  wind.  Approximately  12  hours  after  precipitation  onset,  however  a 
moderate  wind  event  followed  immediately  by  another  of  magnitude  U  ~  5.0m/ s  oc- 
curred. These  increases  in  wind  stress  occurred  during  the  precipitation  event  and  re- 
sulted in  two  increases  of  the  differential  MLD  and  temperature  (points  a  and  b).  Thus, 
the  layer  deepened,  due  to  decreased  We,  causing  the  layer  to  cool. 

Immediately  following  events  a  and  b  the  wind  decreased  to  near  zero  causing 
a  large  decrease  in  MLD.  Here  the  differential  MLD  reached  a  maximum,  as  the  MLD 
in  case  1  shallowed  to  near  the  surface.  Associated  with  this  decrease  in  MLD  in  case 
1,  was  a  negative  differential  temperature.  This  was  the  result  of  the  shallower  MLD  in 
case  1  allowing  for  the  increased  heat  loss  by  the  layer  with  Q0  >  0  (net  heat  loss)  in 
winter  (point  c).  The  differential  MLD  increased  then  decreased,  as  induced  by  the 
onset  and  drop-off  of  a  wind  event  (points  d  and  e).  Again  this  produced  an  associated 
positive  and  negative  differential  temperature,  respectively.  At  day  six  a  strong  storm 
occurred,  (U~  13.0m/s)  (event  f)  followed  immediately  by  a  series  of  larger  fluctuations 
in  wind  speed.  These  storms  brought  the  MLD  of  case  1  to  near  that  of  case  2.  In- 
creased buoyancy,  causing  decreased  entrainment.  caused  the  MLD  in  case  1  to  deepen 
less  than  that  of  case  2  accounting  for  this  significant  differential  temperature  event  at 
day  six  as  case  1.  This  entrained  less  cold  water  than  in  case  2.  After  day  six  a  series 
of  winds  (U  ~  S.Omjs)  continued  to  mix  the  layer,  bringing  the  differential  MLD's  and 
temperatures  closer  to  zero.  Although  the  effect  of  the  precipitation  event  on  the  mixed 
layer  decreased  with  time,  the  layer  temperature  and  depth  were  effected  measureably  for 
this  experiment  period  of  60  days. 

There  was  an  immediate  decrease  in  salinity  as  precipitation  diluted  the  layer, 
and  this  decrease  continued  until  cessation  of  the  precipitation  event.     Mixing  by 
entrainment  decreased  the  differential  salinity  to  a  value  of  -0.4ppm,  after  which  time  it 
rapidly  approached  that  of  case  2,  aided  by  wind  mixing. 
2.     Event  #2  (Jan  1969) 

As  in  Dec  1968,  prior  to  the  precipitation  event  the  MLD  was  relatively  deep 
(h  ~  90m)  and  continued  to  deepen  during  this  winter  regime.    The  values  for  Q0  still 
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indicated    a    net    loss    to    the    atmosphere    and    strong    winds    were    predominate 
(U  ~  5.0m/ m)  . 

After  the  onset  of  the  precipitation  event,  the  results  of  the  event  were  very 
similar  to  those  of  event  #1.  Therefore,  no  specific  figures  were  provided  other  than 
those  found  between  days  30  and  60  in  Fig.  6  on  page  22. 

At  the  beginning  of  the  precipitation,  differential  MLD  generated  a  large  de- 
crease ( ~  —  100.0m)  as  the  MLD  of  case  1  nearly  reached  the  ocean  surface.  This  was 
concurrent  with  a  wind  speed  lull  for  two  days,  enhancing  this  MLD  shallowing.  Asso- 
ciated with  this  MLD  decrease  was  a  decrease  in  the  differential  temperature  as  the 
shallower  MLD  in  case  1  enhanced  the  temperature  decrease  associated  with  a  net  loss 
of  heat  in  winter,  Q0  <  0.  By  day  40  a  series  of  strong  wind  events  (U  ~  l.Smjs)  acted 
to  converge  the  MLD  of  case  1  back  to  near  that  of  case  2.  Likewise,  the  differential 
temperature  approached  zero  in  a  stepwise  fashion  as  the  temperature  in  case  1  ap- 
proached that  of  case  2.  Decreased  entrainment  rate  due  to  increased  buoyancy  did  not 
allow  the  MLD  of  case  1  to  remain  as  deep  as  that  of  case  2.  This  caused  the  temper- 
ature of  case  1  at  day  33  to  exceed  that  of  case  2.  It  did  not,  however,  exceed  the  in- 
crease seen  in  Event  #1  as  the  difference  in  the  MLD's  between  case  1  and  case  2  were 
smaller  in  this  event. 

There  was  an  immediate  decrease  in  salinity  as  precipitation  diluted  the  layer 
and  this  decrease  continued  until  the  end  of  the  rain  event.    Mixing  by  diffusion  and 
advection  decreased  the  differential  salinity  to  a  value  of  -0.3ppm,  after  which  time  it 
rapidly  approached  that  of  case  2. 
3.     Event  #3  (Feb  1969) 

The  MLD  prior  to  the  precipitation  event  was  deep  (h  ~  95m)  and  deepening 
during  the  winter  regime.  The  value  of  Q0  still  indicated  a  net  loss  of  heat  to  the  at- 
mosphere, but  the  magnitudes  of  Q0  and  of  the  wind  speeds  were  smaller  than  for  the 
previous  two  events. 

After  onset  of  the  precipitation  event  in  case  1,  the  results  of  the  event  were  very 
similar  to  those  of  Event  #1  and  #2.  As  in  Event  #2,  no  specific  figures  were  provided 
other  than  those  found  between  days  60  and  90  in  Fig.  6  on  page  22.  As  the  rain  began, 
differential  MLD  generated  a  large  negative  event  (  ~  —  70.0m)  as  the  MLD  in  case  1 
nearly  shallowed  to  the  ocean  surface.  This  was  concurrent  with  a  wind  lull  of  four  days, 
enhancing  this  shallowing.  Associated  with  this  MLD  decrease  was  a  decrease  in  the 
differential  temperature  as  the  shallower  MLD  in  case  1  increased  the  temperature  drop 


28 


associated  with  a  net  surface  heat  loss.  During  the  end  of  the  precipitation  event  there 
was  a  small  peak  in  both  the  differential  MLD's  and  temperatures  associated  with  a  wind 
speed  maximum.  By  day  70  a  series  of  strong  wind  events  (U  ~  5.0m/ s)  caused  the  MLD 
of  case  1  to  converge  to  near  that  of  case  2.  Likewise,  the  differential  temperature  ap- 
proached zero  in  a  stepwise  fashion  as  the  temperature  in  case  1  approached  that  of  case 
2.  Decreased  entrainment  rate  due  to  increased  buoyancy  did  not  allow  the  MLD  of 
case  1  to  remain  as  deep  as  that  of  case  2.  This  caused  temperature  of  case  1  at  day  70 
to  exceed  that  of  case  2.  It  did  not,  however,  exceed  the  increase  seen  in  Events  #1  and 
#2,  as  the  difference  in  the  MLD's  between  case  1  and  case  2  were  smaller  in  this  event. 

There  was  an  immediate  decrease  in  salinity  in  the  first  case,  as  in  the  first  two 
events,  as  the  precipitation    diluted  the  layer.    This  decrease  continued  until  cessation 
of  the  precipitation  event.    Mixing  by  entrainment  decreased  the  differential  salinity  to 
a  value  of  -0.3ppm,  after  which  it  rapidly  approached  that  of  case  2. 
4.     Event  #4  (Mar  1969) 

Examining  the  MLD  model  output  for  the  E-P  =  0  case  (Fig.  6  on  page  22),  it 
appears  that  this  was  the  period  of  deepest  MLD  (h  >  100m).  This  was  also  the  begin- 
ning of  the  period  of  transition  from  the  colder  winter  regime  to  the  warmer  summer 
regime.  Here  the  MLD  began  to  shallow  due  to  decreasing  wind  speed  and  Q0  >  0. 

This  event  provided  different  results  than  Events  #1  to  #3,  and  Fig.  9  on  page 
30  is  presented  for  better  resolution.  As  in  the  previous  events,  there  was  an  initial  de- 
crease in  both  the  differential  MLD  and  in  temperature  (point  a).  This  was  associated 
with  the  precipitation  event  and  a  large  decrease  in  U  from  18.0  to  near  0.0  m's.  During 
this  shallowing  period  there  was  a  slight  wind  increase  causing  U  to  peak  to  5.0  m,  s,  and 
a  deepening  in  both  case  1  and  case  2.  This  was  less  in  case  1  than  in  case  2.  This  dif- 
ferential in  MLD  generated  a  positive  differential  temperature  (point  b).  The  positive 
differential  temperature  immediately  decreased  as  the  wind  dropped  to  zero,  and  the 
MLD  became  very  small,  order  1  m  (point  c).  Net  heat  loss  cooled  the  upper  layer  in 
case  1  more  than  in  case  2,  generating  a  severe  negative  spike  in  the  differential  tem- 
perature ( 0.10°C).  A  series  of  mild  wind  speeds  until  day  98  yielded  various  negative 

and  positive  peaks  in  the  differential  MLD's  and  temperatures. 

At  day  100  a  strong  wind  (U  ~  8.0m/.s)  returned  the  MLD  in  case  1  equal  to  that 
of  case  2,  making  the  differential  MLD  zero  (point  d).  This  mixing  decreased  the  dif- 
ferential temperature  to  near  zero,  but  it  remained  slightly  colder  as  the  entrained  water 
from  below  the  mixed  layer  remained  warmer  than  water  of  the  previous  shallow  layer. 
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Strong  winds  at  days  113,  130,  and  140  continued  entrainment.  thus,  diffusing  the  colder 
layer  in  case  1  with  the  warmer  water  below,  and  bringing  the  differential  temperature 
closer  to  zero  (points  e  to  g). 

Near  the  end  o^this  60-day  experiment  there  was  a  very  strong  mixing  event  at 
day  147,  induced  by  a  U  peak  of  15.0  m/s  with  a  duration  of  3  days.  Almost  60  days 
after  the  single  precipitation  event,  there  was  still  an  effect  on  the  mixed  layer.  Exam- 
ining the  differential  MLD  and  temperature,  it  appeared  that  the  MLD  in  case  1  did  not 
deepen  as  much  as  in  case  2.  This  deeper  MLD  caused  the  temperature  in  case  1  to  begin 
to  exceed  that  of  case  2  (point  h). 

Following  the  cessation  of  precipitation,  in  the  first  case,  diluted  salinity  was 
mixed  by  entrainment,  decreasing  the  differential  salinity  to  a  value  of  -0.4ppm.  After 
this  it  rapidly  approached  that  of  case  2  by  wind  mixing. 
5.     Event  #5  (Apr  1969) 

Near  the  end  of  the  March  study  period  (Event  #4)  the  winds  decreased  and  the 
MLD  shallowed  (Fig.  6  on  page  22).  After  the  wind  speed  lull  from  day  123  to  day  127, 
there  was  a  steady  rise  in  wind  to  the  peak  value  at  day  147  of  U~  14.0ml s  .  This 
gradual  increase  provided  the  TKE  to  mix  the  layer  to  to  the  deepest  point  for  the  year 
(h  ~  125.0m).  This  was  followed  by  a  period  of  shallowing,  induced  by  an  abrupt  fall  in 
wind  speed  (Fig.  10  on  page  32).  The  value  of  Q0  was  positive  and  increasing,  net  heat 
flux.  This  served  to  provided  a  net  heat  gain  in  this  summer  regime,  as  opposed  to  the 
heat  loss  that  occurred  with  the  winter  regime. 

Fig.  10  on  page  32  demonstrated  the  initial  decrease  of  MLD  in  case  1,  to  be 
greater  than  that  of  case  2  (point  a).  Due  to  the  low  wind  speeds  coincident  with  the 
initiation  of  the  precipitation  event,  the  mixed  layer  of  case  2  was  also  relatively  shallow 
(h  ~~  21.0m)  .  Thus  the  increase  in  the  differential  temperature  was  less  than  that  of  one 
day  later  (at  point  b)  following  a  wind-driven  deepening  event.  This  produced  a  large 
decrease  in  the  differential  MLD.  A  series  of  shallowing  and  deepening  events  followed 
and  were  associated  with  differential  temperature  adjustments.  This  brought  the  differ- 
ential MLD  and  the  differential  temperature  back  toward  zero. 

The  largest  decrease  in  differential  salinity  in  Fig.  6  on  page  22  occurred  in  this 
event.  In  this  event  the  precipitation  event  was  initiated  after  a  long  period  of  low 
winds,  providing  a  shallow  MLD  prior  to  onset  of  the  precipitation  event.  This  may 
have  been  an  unrealistic  example.  Given  that  increased  wind  speeds  generally  ac- 
company large  precipitation  events  like  that  of  day  124,  this  hypothetical  event  would 
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have  been  unlikely  to  occur,  but  it  permits  an  interesting  examination  of  model  dynam- 
ics. 

6.  Event  #6  (May  1969) 

Between  day  150  and  day  180  the  MLD  in  case  2  was  similar  to  that  of  Event 
#5,  but  more  shallow.  This  shallowing  was  the  result  of  decreased  wind  forcing  as  seen 
in  Fig.  6  on  page  22.  Net  heat  gain  to  the  ocean  surface  continued  to  increase  as  sum- 
mer approached,  and  some  of  the  largest  net  downward  heat  fluxes  occurred  during  this 
period  (&  miX~  700  W). 

The  dynamics  of  this  event  resembled  those  of  Event  #5.  No  60-day  period  fig- 
ures were  provided.  Fig.  6  on  page  22  provides  adequate  resolution.  As  before,  there 
was  an  initial  shallowing  producing  a  positive  differential  temperature  associated  with 
the  shallower  MLD  and  the  summer  heating  regime.  Until  day  165,  the  differential 
MLD  and  the  differential  temperature  resembled  those  of  the  previous  events.  After  day 
165,  however,  the  resulting  differential  temperature  profile  differed  from  those  of  the 
previous  events.  Here,  the  differential  temperature  reached  a  post-precipitation  mini- 
mum value  of  0.05°C  and  continued  to  increase  slowly  through  the  rest  of  the  study 
period. 

This  steady  increase  in  differential  temperature,  when  compared  with  the  tem- 
perature plot  of  the  mixed  layer,  indictated  that  the  mixed  layer  was  warming  faster  for 
case  1  than  case  2.  The  MLD's  were  equal  in  depth  during  this  period,  so  that  this 
condition  could  not  have  been  due  to  case  1  being  shallower  than  case  2.  The  possible 
causes  were  decreased  \Ve,  due  to  decreased  salinity,  and,' or  increased  AT.  An  increased 
AT  could  be  caused  by  the  mixed  layer  in  case  1  being  warmer  than  that  of  case  2,  and 
the  temperature  of  the  water  mass  below  being  the  same  for  both  cases.  If  so,  in  the 
summer  regime  this  steady  increase  in  the  differential  temperature  would  appear  to  have 
been  "self-perpetuating." 

There  was  an  immediate  decrease  in  salinity  in  the  first  case,  as  in  the  other 
events,  as  the  precipitation  diluted  the  layer.  This  decrease  continued  until  cessation 
of  the  precipitation  event  at  day  156.  Mixing  by  entrainment  decreased  the  differential 
salinity  to  a  value  of  -0.3ppm,  after  which  time  it  rapidly  approached  that  of  case  2. 

7.  Event  #7  (Jun  1969) 

Though  the  MLD  for  case  2  began  relatively  deep  (h  ~  70m)  it  decreased  to  the 
shallowest  mean  period  value  (h  ~  20m).  This  was  caused  by  an  initial  wind  event  at  day 
181  {U~6mjs)  ,  and  followed    by  a  relative  lull  from  day  182  to  205.    Another    wind 
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Fig. 


11.      Event  #7  (Jim  1968):     (a)  differential  MLD,  (b)  differential  temperature, 
(c)  differential  salinity,  and  (d)  wind  speed 
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increase  occurred,  similar  in  magnitude  to  the  first  at  day  205  (Fig.  1 1  on  page  34).  The 
net  heat  flux  remained  strongly  downward,  but  not  as  strong  as  in  the  previous  study 
(Fig.  7  on  page  23). 

Initially,  the  MLD  for  case  2  was  shallow  (h2  ~  11m)  .  At  the  beginning  of 
precipitation,  the  MLD  of  case  1  shallowed  to  only  3m,  yet  there  remained  a  small  dif- 
ferential MLD,  producing  a  small  differential  temperature  <0.1°C  (Fig.  11  on  page 
34,point  a).  This  difference  decreased  with  continuing  low  level  winds  which  acted  to 
deepen  the  MLD  in  both  cases  and  yielded  a  decrease  in  the  differential  temperature  to 
near  zero  (point  b). 

The  most  significant  deepening  event  of  this  period  then  occurred  (point  c), 
which  produced  the  largest  differential  temperature  of  the  period  (  ~  0.5° Q  .  The  mixed 
layer  of  the  rain  case  then  behaved  as  expected,  stepping  differential  MLD  and  differ- 
ential temperature  back  toward  zero.  The  differential  temperature  for  this  period  was 
much  larger  than  for  the  other  periods  studied.  This  appeared  to  be  due  to  the  very 
shallow  MLD's  achieved  by  case  1,  enhancing  the  effect  of  the  downward  summer  heat 
flux. 

The  shallow  MLD  at  the  onset  of  precipitation  did  yield  a  stronger  than  normal 
differential  salinity  (  ~  —  Q.QIQppm).    The  differential  salinity  acted  as  in  the  previous 
events,  returning  to  an  equilibrium  value  of  ~  —  0.02ppm  . 
8.     Event  #8  (Jul  1969) 

This  event  coincided  with  a  shallow  (h  ~  30m)  trend  in  the  MLD  for  case  2. 
This  was  associated  with  low  wind  speeds  (U  <  5ml s)  (Fig.  7  on  page  23  and  Fig.  12  on 
page  36).   Net  heat  flux  appeared  to  be  steady  (<90  ~  180  W)  . 

There  was  an  initial  positive  increase  in  differential  temperature  which  occurred 
in  a  stepwise  manner,  commensurate  with  the  differential  MLD  (Fig.  12  on  page  36). 
This  occurred  until  wind  speed  decreased  to  zero  at  day  216  (point  a).  From  day  216 
to  226  (points  a  to  b)  the  differential  temperature  remained  nearly  constant  at 
~  0.065° C,  with  only  a  slight  increase  in  time.  The  MLD  in  case  1  equaled  that  of  case 
2  (and  the  differential  MLD  ~  0)  in  this  period.  This  condition  was  similar  to  the  slight 
increase  which  was  notel  in  for  event  #6,  and  the  physical  explanation  is  the  same. 

After  day  226  (point  b)  there  was  a  large  increase  in  the  differential  MLD  and 
differential  temperature  due  to  a  U  spike  of  ~  4.3m/ s.  This  acted  to  cool  the  mixed  layer 
of  case  2  more  than  that  of  case  1  and  yielded  a  positive  increase  in  the  differential 
temperature  (  ~  0.1° C).    The  differential  temperature  then  fluctuated  around  0.22°C  for 
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the  remainder  of  this  60-day  period  as  mild  winds  allowed  shallowing  of  the  mixed  layers 
in  cases  1  and  2  with  intermittent  weak  spikes  in  wind  speed  producing  some  oscillation 
around  this  mean. 

There  was  an  immediate  decrease  in  salinity  in  the  first  case,  as  in  the  other 
events,  as  the  precipitation  diluted  the  layer  and  this  decrease  continued  until  cessation 
of  the  precipitation  event.  iMixing  by  entrainment  decreased  the  differential  salinity  to 
a  value  of  -0.3ppm,  after  which  time  it  rapidly  approached  that  of  case  2. 

9.  Event  #9  (Aug  1969) 

With  the  exception  of  the  initial  small  increase  in  differential  temperature,  this 
event  produced  results  very  much  like  that  of  Event  #8.  This  small  fluctuation  was  ac- 
counted for  by  the  weak  wind  mixing  (U  ~  4ml s)  occurring  upon  onset  of  precipitation. 

10.  Event  #10  (Sep  1969) 

During  this  study  period  the  MLD  in  case  2  was  relatively  shallow,  and  begin- 
ning to  deepen,  as  there  was  a  seasonal  shift  from  summer  to  winter  (Fig.  7  on  page 
23).  This  was  associated  with  steadily  increasing  wind  speeds  (Fig.  13  on  page  38).  As 
the  seasonal  regime  began  to  shift  to  a  winter  regime,  the  Q0  shifted  from  one  of  net 
downward  heating  toward  one  of  cooling,  although,  during  this  period  net  heating  re- 
mained the  slightly  positive. 

As  with  the  results  of  Event  #9,  these  experimental  results  resembled  those  of 
Event  #8  without  the  early  fluctuating.  In  this  event,  prior  to  precipitation  onset,  the 
MLD  in  both  cases  had  shallowed  to  near  surface  due  to  a  wind  speed  lull  from  day  273 
to  275  (Fig.  13  on  page  38).  A  slight  fluctuation  at  day  274  (point  a)  generated  a  small 
downward  event  in  differential  MLD  (~  —  \m)  which  promoted  a  slight  upward  differ- 
ential temperature.  With  a  follow  on  U  event  at  day  275  (point  b)  of  ~  4ml s,  the  largest 
of  all  positive  differential  temperature  in  the  twelve  events  was  generated.  This  occurred 
along  with  a  differential  MLD  of  ~  -  25m.  Wind  events  followed  which  increasingly 
deepened  the  MLD  of  case  1  and  brought  the  differential  temperature  back  to  near  zero, 
and  constantly  decreasing. 

This  event  produced  the  largest  differential  salinity  of  all,  concurrent  with  this 
anomalous  differential  temperature  fluctuation.  Again,  the  very  shallow  MLD  estab- 
lished before  precipitation  onset  appeared  to  have  contributed  to  this  effect,  and  equi- 
librium was  quickly  reestablished  at  differential  salinity  ~  —  O.Olppm. 


37 


a .  0 

i  8 


o 

I 


£<= 


b. 


c 

o 


IT) 
O 


^d. 


3 


270 


C. 


^  o 


<0 


to 


d. 


CM 


o 


CD 

z:  o 


1 — 

270     280 

b 


280 


270     280 


290 


290 


300 


300 


310 


310 


290      300      310 

TIMEtDYJ 


270     280 


290 


300 


310 


320 


320 


320 


330 


330 


330 


b 


320 


330 


Fig.  13.      Event  #10  (Sep  1969):    (a)  differential  MLD,  (b)  dillcrential  temperature, 
(c)  dillcrential  salinity,  and  (d)  wind  speed 


38 


11.  Event  #11  (Oct  1969) 

In  this  period  the  MLD  of  case  2  was  becoming  deep  again  (h  ^  60m)  and 
continued  to  deepen  in  this  early  winter  regime  (Fig.  7  on  page  23).  Increasing  wind 
speeds  and  near  zero  net  heat  values  accounted  for  this  elTect. 

This  return  to  the  winter  regime  bidded  a  return  to  the  differential  MLD  and 
differential  temperature  profiles  seen  in  the  first  three  events.  Initial  cooling  of  the 
shallower  MLD  of  case  1  generated  a  negative  spike  in  the  differential  temperature 
which  moved  positive  as  the  decreased  entrainment  rate  of  case  2  inhibited  entrainment 
of  the  layer.  Further  mixing  acted  to  equalize  MLD's  and  return  the  differential  MLD 
and  differential  temperature  toward  zero.  After  the  differential  temperature  flux  moved 
to  a  positive  value  it  remained  at  a  larger  level  than  observed  in  the  first  three  events, 
as  it  did  in  the  summer  regimes.  This  indicated  that  the  mixed  layer  was  in  transition 
from  the  summer  to  winter  regimes  and  MLD's  were  shallower  than  those  of  the  first 
cases. 

There  was  an  immediate  decrease  in  salinity  in  the  first  case,  as  in  the  other  ex- 
periments, as  the  precipitation  diluted  the  layer  and  this  decrease  continued  until  ces- 
sation of  the  precipitation  event.  Mixing  by  entrainment  decreased  the  differential 
salinity  to  a  value  of  -0.2ppm,  after  which  time  it  rapidly  approached  that  of  case  2. 

12.  Event  #12  (Nov  1969) 

This  event  was  very  similar  in  forcing  conditions  and  resulting  output  to  that 
of  Event  #11.  In  this  event  the  MLD  in  case  2  was  deeper  than  that  of  Event  #11  as 
early  winter  deepening  continued,  and  mean  wind  speeds  appeared  to  be  stronger  with 
larger  peak  wind  events  (Fig.  7  on  page  23).  This  deeper  MLD  yielded  a  larger  differ- 
ential MLD  after  precipitation  onset,  as  well  as  a  lesser  positive  peak  of  differential 
temperature.   This  also  yielded  a  similar  differential  salinity  perturbation. 

C.     EXPERIMENT  SUMMARY 

The  results  of  these  events  indicate  that  the  effects  of  discrete  precipitation  events 
varied  greatly  as  a  function  of  low  frequency  seasonal  heat  flux  forcing  and  higher  fre- 
quency wind  forcing.  These  two  forcing  boundary  conditions  strongly  influenced  the 
mixed  layer  thermal  characteristics,  as  well  as  the  salinity. 

The  average  deeper,  colder  mixed  layer  in  winter  was  significantly  shallowed  by  a 
single  precipitation  event.  This  shallowing  allowed  for  the  increased  loss  of  the  thermal 
energy  of  the  shallower  mixed  layer  to  the  colder  atmosphere,  supporting  a  temperature 
decrease  in  the  mixed  layer.   Wind  events  that  followed  precipitation  deepened  the  mixed 
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layer,  but  it  did  not  deepen  as  much,  nor  as  quickly  as  in  the  case  without  precipitation. 
Follow-on  wind  events  deepened  the  mixed  layer  until  entrainment  returned  the  MLD 
to  that  of  the  case  without  precipitation.  Entrainment  had  a  similar  effect  on  the  dif- 
ferential mixed  layer  temperature. 

The  effects  of  the  precipitation  event  in  winter  were  modulated  by  the  discrete 
deepening  and  shallowing  events  that  occurred  just  prior  to  precipitation  onset.  A 
deeper  MLD  increased  the  magnitudes  of  the  differential  values  that  occurred  as  a  result 
of  precipitation  shallowing.  A  shallower  MLD  operated  in  the  opposite  sense  to  de- 
crease these  values.  Deepening  and  shallowing  events  occurring  concurrent  with  pre- 
cipitation had  a  similar  effect,  acting  to  enhance  or  decreasing  differential  values. 

In  summer  the  average  shallower,  warmer  mixed  layer  was  less  significantly  shal- 
lowed by  these  precipitation  events.  Given  that  the  mixed  layer  was  already  relatively 
shallow,  the  decrease  in  MLD  was  not  as  significant  as  in  winter.  Significant  thermal 
stratification  in  summer  yielded  a  large  increase  in  mixed  layer  temperature  with  precip- 
itation induced  shallowing,  as  the  thermal  energy  input  at  the  top  of  the  mixed  layer  was 
concentrated  into  a  shallower  mixed  layer.  This  increased  mixed  layer  temperature 
produced  a  more  stable  mixed  layer  as  the  temperature  jump  at  the  base  of  the  mixed 
layer  increases.  The  increase  in  temperature  jump  tended  to  further  inhibit  entrainment 
and  slowed  deepening  of  the  mixed  layer.  Deepening  slowed  enough  that  the  differential 
temperature  actually  increased  with  time. 

The  discrete  precipitation  events  could  have  had  a  long  term  effect  on  the  mixed 
layer.  This  was  especially  true  in  summer  when  MLD's  were  shallower.  In  cases  where 
strong  winds  generated  strong  mixing  events,  the  mixed  layer  continued  to  demonstrate 
increased  buoyancy,  as  long  as  55  days  after  the  precipitation  event,  as  the  differential 
MLD  became  negative  and  differential  temperature  yielded  a  positive  fluctuation. 
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V.     PROCEDURES,  RESULTS,  AND  ANALYSIS 

A.     PROCEDURES 

The  400-day  mixed  layer  simulation,  as  described  previously,  was  conducted  to  ex- 
amine the  accumulative  effect  of  a  time  series  of  rainfall  events  in  explaining  salinity  and 
temperature  changes  with  time  and  depth.  Discrete  daily  precipitation  data  was  not 
available  for  this  study.  Only  monthly  totals  were  available.  As  a  result  a  method  was 
devised  to  simulate  daily  precipitation  values  using  the  available  oceanographic  and 
meteorological  observations.  With  the  lack  of  surface  pressure  and  upper  atmospheric 
data,  precipitation  prediction  became  difficult,  thus  prediction  of  localized  atmospheric 
convective  activity  was  difficult.  Utilizing  wind  and  cloud  information,  however,  an  at- 
tempt has  been  made  to  predict  precipitation  associated  with  frontal  activity.  This 
technique  was  based  on  meteorological  principles  developed  by  Petersen  (1956). 

Specifically,  this  procedure  utilized  wind  speed  and  direction,  and  cloud  coverage 
data  (Fig.  14  on  page  42).  With  cloud  cover  data  available  in  octaves,  precipitation  was 
allowed  for  cases  with  cloud  coverage  equal  to  or  greater  than  7/8.  Considering  pre- 
frontal and  frontal  rain,  precipitation  was  allowed  to  occur  for  cases  with  wind  directions 
from  000°  (northerly  Wi:ids)  clockwise  to  220°  (southwesterly  winds),  and  wind  speeds 
greater  than  10  kts.  This  accounted  for  the  southerly  to  southwesterly  winds  that  pro- 
ceed a  front  and  the  "confused"  northerly  to  southeasterly  winds  that  accompany  frontal 
passage.  For  post-frontal  rain,  precipitation  was  allowed  to  occur  for  cases  with  wind 
direction  220°  (southwesterly  winds)  clockwise  to  000°  (northerly  winds),  and  for 
stronger  wind  speeds,  greater  than  25  kts.  This  accounted  for  the  strong  southwesterly 
to  northwesterly  winds  that  accompany  a  post-frontal  regime  and  ensured  consideration 
of  all  possible  wind  directions.  Additionally,  the  post-frontal  case  required  that  the  sea 
surface  temperature  be  greater  than  the  air  temperature,  as  the  cold  frontal  air  passes 
over  warmer  water.  Precipitation  amount  was  adjusted  to  yield  the  400-day  precipitation 
total  of  the  observed  monthly  data.  The  resulting  rain  with  time  is  contained  in  Fig.  15 
on  page  43. 

This  "complex"  alogarithm  for  precipitation  forcing  scheme  was  applied  to  the 
OPBL  model  at  a  range  of  possible  amounts  that  may  represent  normal  interannual 
variability  in  the  North  Pacific.  A  "medium"  precipitation  level  was  chosen  by  adjusting 
total  400-day  output  to  equal  that  of  the   summation  of  the  13  month  period  from 
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Fig.   14. 


Front  and  associated  wind  directions  and  isotherms 
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Fig.  15.  Complex  Rain  Model  Output  and  Winds  for  400  day  period:  (a)  rain 
model  output  for  400  days  from  23  Nov  1968  adjusted  to  equal  actual 
monthly  totals,  and  (b)  wind  for  same  time  period. 
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December  1968  to  January  1970  (Table  3  on  page  45).  The  "low"  level  was  chosen  to 
equal  two-third's  of  the  total  400-day  precipitation  total  of  the  medium  level  case.  The 
"high"  level  case  was  equal  to  one  and  one-half  times  the  medium  case  and  the  "very 
high"  case  was  made  equal  to  three  times  the  medium  case  (Fig.  16  on  page  46).  The 
same  wind  stress  and  heat  flux  atmospheric  forcing  values  and  initial  conditions  de- 
scribed in  chapter  II  were  used  for  all  four  model  runs  as  in  chapter  IV. 

Table  2.     PRECIPITATION  TOTALS  FOR  THE  400-DAY  STUDY  PERIOD 


Precipitation  Forcing  Totals  for  400-day  Period 

qualitative  designation 

quantitative  400-day  total 

low 

32  cm 

medium 

48  cm 

high 

96  cm 

very  high 

160  cm 

These  four  complex  cases  were  compared  with  the  effects  of  "simple"  constant 
evaporation  and  precipitation.  Here  the  four  "simple"  precipitation  and  evaporation 
rates  were  specified  to  give  constant  evaporation  and  precipitation  with  the  same  total 
precipitation  and  evaporation  as  those  of  the  four  complex  cases.  These  simple  precip- 
itation cases  were  applied  to  the  model  so  that  the  effects  of  discrete  precipitation  events 
could  be  compared  to  those  of  constant  E-P.  The  resulting  MLD,  temperature,  and 
salinity  time  series  were  differenced  for  comparison.  The  output  mixed  layer  values  o[ 
MLD,  mixed  layer  temperature,  and  mixed  layer  salinity  from  complex  and  simple  forc- 
ing were  compared  to  the  CTD  observations. 

Advection  was  maintained  in  the  model  at  0.0  ppm/day.  Considering  the  placement 
of  advection  in  the  salinity  equation: 

dS       S(E-P)      ASWe 


tit 


It  should  be  noted  that  an  increase  in  advection  would  provide  the  same  tendency  for 
mixed  layer  salinity  as  would  decreasing  precipitation. 
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Table  3.     MONTHLY  PRECIPITATION  AT  OCEAN  STATION  "P"  FOR  400 
DAY  PERIOD 


Monthly  Precipitation  for  400-day  Period 

month/year 

precipitation 

Dec:  1968 

1.5cm 

Jan  1969 

0.6cm 

Feb  1969 

2.9cm 

Mar  1969 

3.3cm 

Apr  1969 

4.2cm 

May  1969 

3.1cm 

Jun  1969 

2.1cm 

Jul  1969 

1.3cm 

Aug  1969 

1.3cm 

Sep  1969 

3.1cm 

Oct  1969 

5.5cm 

Nov,  1969 

8.0cm 

Dec  1969 

6.3cm 

B.     RESULTS  AND  ANALYSIS 

1.     Complex  Precipitation  and  Evaporation  Forcing 

The  response  of  the  OPBL  model  depended  upon  the  net  precipitation  for  the 
cases  of  complex  precipitation.  The  400-day  total  precipitation  was  increased  from  low 
to  very  high  as  indicated  in  Table  2  on  page  44.  An  increase  in  total  precipitation  caused 
a  general  decrease  in  MLD's,  as  a  result  of  increased  buoyancy  and  decreased 
entrainment  (Fig.  17  on  page  47).  As  the  precipitation  level  for  complex  precipitation 
forcing  was  increased  from  low  to  very  high,  it  was  larger  in  magnitude  in  winter  than 
in  summer,  though  the  relative  increase  in  summer  was  proportional  to  that  of  winter 
(Fig.  16  on  page  46).  This  was  due  to  the  small  amount  of  precipitation  in  summer 
compared  to  winter.  This  increase  in  precipitation  in  winter  with  increasing  total  pre- 
cipitation increased  buoyancy  more  in  winter  than  in  summer,  yielding  shallower  MLD's 
for  larger  complex  precipitation  amounts.  Additionally,  the  largest  magnitude  of  differ- 
ential MLD's  occurred  in  winter  (Fig.  17  on  page  47  and  Fig.  18  on  page  48). 

Associated  with  the  differential  MLD  flux  that  accompanied  increased  complex 
precipitation  levels  was  a  change  in  mixed  layer  temperature  and  salinity.    Differential 
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Fig.  16.  Simple  and  Complex  Precipitation  Forcing:  (a)  low  level  precipitation, 
(b)  medium  level  precipitation,  (c)  high  level  precipitation,  and  (d)  very 
high  precipitation. 
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Fig.  17.  Complex  Precipitation  Forcing  Predicted  MLD's:  (a)  solid  line  represents 
very  high  precipitation  forced  case  and  (b)  dotted  line  represents  low  pre- 
cipitation forced  case. 
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Fig.  18.  Subtracted  Values  of  Complex  Precipitation  Forcing  Output:  All  values 
arc  the  result  of  subtracting  Very  High  precipitation  forcing  output  from 
Low  precipitation  forcing,  (a)  differential  MLD,  (b)  differential  temper- 
ature, and  (c)  differential  salinity. 
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Fig.  19.  Complex  Precipitation  Forcing  Predicted  Temperatures:  (a)  dotted  line 
was  low  precipitation  forced,  (b)  dashed  line  was  medium  precipitation 
forced,  (c)  chained-dashed  was  high  precipitation  forced,  and  (d)  solid  line 
was  very  high  precipitation  forced. 
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mixed  layer  temperature  decreased  only  slightly  in  winter  when  the  MLD  was  deep. 
Negative  differential  temperatures  were  induced  by  the  winter  heat  loss  of  the  shallower 
mixed  layer.  As  the  summer  approached,  and  MLD  shallowed  due  to  increased  Q0  and 
decreased  winds,  the  differential  MLD,  had  a  more  significant  effect  on  differential 
temperature.  Thus,  the  decreased  MLD  in  the  shallowing  season  concentrated  more  of 
the  incident  solar  energy  in  a  shallower  layer  than  it  would  have  without  increased  pre- 
cipitation.  This  produced  a  large  positive  differential  temperature. 

Examination  of  the  differential  salinity  showed  a  continual  decrease  in  salinity 
as  very  high  complex  precipitation  forcing  inputed  more  fresh  water  into  the  mixed  layer. 
Deepening  events  which  yielded  differential  MLD's  produced  differing  entrainment 
events.  Deeper  MLD's,  characteristic  of  lower  complex  precipitation  forced  cases,  en- 
trained water  of  higher  salinity  than  in  the  higher  complex  case.  Larger  precipitation 
events  input  freshwater  into  the  surface  layer  inhibiting  entrainment,  yielding  a  shallower 
MLD.  This  reduction  in  MLD  further  concentrated  the  fresh  water,  reducing  surface 
S  (Fig.  20  on  page  51). 

From  CTD  observations,  actual  values  of  MLD.  temperature,  and  salinity  were 
compared  to  those  values  predicted  by  the  OPBL  model  with  complex  precipitation 
forcing.  Verification  was  hindered  by  data  voids,  especially  during  the  first  95  days  of 
the  study.  Nevertheless,  it  may  be  concluded  that  the  MLD's  of  the  CTD  data  were  best 
approximated  by  the  model  using  very  high  complex  precipitation.  The  envelopes  of  the 
fluctuating  MLD  of  all  four  complex  precipitation  forcing  cases  appeared  to  be  deeper 
than  those  of  the  CTD  casts.  With  the  small  number  of  CTD's  over  any  synoptic  pe- 
riod, an  hour-by-hour  comparison  was  not  possible,  although  MLD  was  predicted  to 
fluctuate  diurnally  by  as  much  as  50  m  day.  As  the  level  of  precipitation  was  increased, 
the  resulting  MLD  became  shallower.  Thus,  the  case  with  the  largest  precipitation 
amount  appeared  to  most  closely  approximate  the  CTD  envelope  for  MLD  (Fig.  21  on 
page  52). 

The  mixed  layer  temperatures  compared  favorably  with  observations  in  all  four 
complex  cases.  As  the  study  period  progressed,  the  temperatures  of  the  actual  data  in- 
creased more  than  those  from  complex  precipitation  forcing  commensurate  with  the  ac- 
tual MLD  being  shallower.  Thus,  the  case  with  the  highest  complex  precipitation 
produced  the  greatest  shallowing,  yielding  the  highest  summer  temperature.  The  peak 
temperature  from  the  very  high  complex  precipitation  case,  however,  only  slightly  ex- 
ceeded that  of  the  low  precipitation  case  (Fig.  22  on  page  53). 
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Fig.  20.  Complex  Precipitation  Forcing  Predicted  Salinities:  (a)  dotted  line  was 
low  precipitation  forced,  (b)  dashed  line  was  medium  precipitation  forced, 
(c)  chained-dashed  was  high  precipitation  forced,  and  (d)  solid  line  was 
very  high  precipitation  forced. 
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Fig.  21.      Complex  Precipitation  Forcing  Predicted  MLD  and  CTD  MLD:     (a)  the 

solid  line  was  very  high  precipitation  forced  and  (b)  the  dotted  line  was 
low  precipitation  forced. 


S2 


Fig.  22.       Complex   Precipitation  Forcing   Predicted  and  CTD  Temperature:     (a) 

dotted  line  was  low  precipitation  forced,  (b)  dashed  line  was  medium  pre- 
cipitation forced,  (c)  chained-dashed  was  high  precipitation  forced,  and  (d) 
solid  line  was  very  high  precipitation  forced. 
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Fig.  23.  Complex  Precipitation  Forcing  Predicted  and  CTD  Salinity:  (a)  dotted 
line  was  low  precipitation  forced,  (b)  dashed  line  was  medium  precipi- 
tation forced,  (c)  chained-dashcd  was  high  precipitation  forced,  and  (d) 
solid  line  was  very  high  precipitation  forced. 
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The  model-predicted  salinity  values  yielded  less  certain  results.  Actual  CTD 
salinity  data  behaved  in  a  cyclic  manner  during  the  study  period.  Here  the  initial  salinity 
of  32.5  ppm  immediately  increased  with  increasing  MLD.  As  the  observed  MLD  shal- 
lowed, the  salinity  decreased,  until  the  minimum  was  achieved  in  summer.  This  was 
possibly  associated  with  shallowing  due  to  decreased  wind  and  precipitation.  As  MLD 
deepened  again  around  day  330,  salinity  increased  (Fig.  23  on  page  54). 

Overall,  this  observed  cyclic  behavior  was  best  approximated  by  modeling  with 
high  complex  precipitation  forcing.  The  two  cases  with  precipitation  less  than  the  high 
case  yielded  salinity  values  that  consistently  exceeded  those  of  the  actual  data,  and  were 
non-cyclic.  During  winter  deepening  for  these  cases,  the  salinity  increased  dramatically 
from  day  0  to  day  100.  During  spring  shallowing  and  into  summer,  the  salinity  contin- 
ued to  increase  gradually,  and  did  not  exhibited  a  cyclical  decrease.  For  the  very  high 
complex  case,  the  salinity  exhibits  a  decreasing  trend  almost  immediately.  Here  the 
buoyant  damping  of  the  turbulence  by  the  strong  precipitation  generally  countered  the 
shear  production  of  turbulence  attributable  to  the  wind  stress. 

Though  the  high  complex  precipitation  case  did  agree  with  the  cyclic  behavior 
of  the  CTD  salinity  data,  there  were  salinity  observations  of  significant  disagreement 
with  the  model  prediction.  As  stated  previously,  peak  salinity  in  winter  was  less  than  the 
actual  value.  In  summer  the  salinity  minimum  was  larger  than  the  observed  minimum 
by  ~  0.7  ppm. 

2.     Constant  Precipitation  and  Evaporation  Forcing 

The  results  obtained  during  this  study  were  similar  to  those  for  the  complex 
precipitation  forcing  study.  Upon  comparing  the  various  magnitudes  of  simple  precipi- 
tation forcing,  it  was  expected  that  increasing  levels  of  constant  E-P  would  have  a 
magnifying  effect  on  the  mixed  layer  parameters.  For  simplicity,  high  and  low  precipi- 
tation forcing  cases  were  compared  here.  The  increased  level  of  precipitation  in  the  high 
case  yielded  a  shallower  overall  MLD  than  the  low  case  (Fig.  24  on  page  56).  This 
shallowing  effect  was  stronger  in  winter  than  in  summer,  as  MLD's  tended  to  be  deeper 
in  winter.  The  peak  predicted  in  the  shallower  period  in  the  low  case  (Fig.  24  on  page 
56  point  a)  did  not  appear  in  the  high  case. 

It  appeared  that  winter  temperatures  were  well  verified  for  all  simple  and  com- 
plex forced  cases  by  the  actual  mixed  layer  temperatures,  as  the  increased  precipitation 
yielded  only  slightly  more  shallowing.     In  summer,  however,  differential  temperatures 
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Fig.  24.  Simple  Precipitation  Forcing  Predicted  MLD'S:  (a)  the  solid  line  was 
very  high  precipitation  forced  and  (b)  the  dotted  line  was  low  precipitation 
forced. 
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Fig.  25.  Simple  Precipitation  Forcing  Predicted  Temperature:  (a)  dotted  line  was 
low  precipitation  forced,  (b)  dashed  line  was  medium  precipitation  forced, 
(c)  chained-dashed  was  high  precipitation  forced,  and  (d)  solid  line  was 
very  high  precipitation  forced. 
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increased  as  the  high  case  absorbed  more  solar  heat  in  a  shallower  layer  than  the  low 
case. 

As  with  the  complex  precipitation  cases,  the  constant  low  and  medium  precipi- 
tation predictions  yielded  a  proportional  decrease  in  salinity.  The  high  precipitation  case 
proved  to  be  the  most  cyclic  in  time  over  the  year.  The  very  high  case  yielded  an  almost 
monotonic  decrease  in  salinity  as  a  result  of  the  large  constant  surface  fresh  water  influx 
(Fig.  26  on  page  59). 

From  comparing  simple  to  complex  precipitation  forcing,  it  was  noted  that  the 
complex  case  produced  relatively  more  precipitation  in  winter  than  in  summer,  and  there 
was  greater  evaporation  in  summer  than  in  winter.  The  constant  E-P  forcing  effect  of 
the  simple  case  thus  had  a  larger  precipitation  rate  during  periods  of  no  rain  for  the 
complex  case  (in  winter  when  low  winds  occur  and  during  the  summer).  In  winter  dif- 
ferential MLD  fluctuated  both  positively  and  negatively  about  zero.  It  appeared  to 
fluctuate  negatively  during  a  strong  complex  precipitation  event  yielding  shallowing.  It 
fluctuated  positively  during  non-rain  periods  for  the  complex  case,  producing  relative 
shallowing  in  the  simple  case.  In  winter  these  differential  MLD  values  exceeded  those 
of  summer. 

Differential  temperature  in  winter  showed  approximately  equal  temperatures  for 
the  four  rain  levels  in  both  the  simple  and  complex  cases.  Additionally,  any  fluctuation 
in  MLD  had  a  relatively  small  effect  on  the  differential  temperature  in  the  winter  season. 
In  summer  a  negative  differential  temperature  was  induced  by  precipitation  being  greater 
in  the  simple  case  than  in  the  complex  case. 

With  respect  to  salinity,  differential  salinity  tended  to  be  positive  in  winter. 
Here,  the  salinity  of  the  complex  case  was  usually  less  than  that  of  the  simple  case  due 
to  the  greater  amount  of  precipitation  in  the  complex  case.  It  did  fluctuate  during  this 
period,  as  periods  of  low  winds  yielded  a  decreased  complex  precipitation  allowing  sim- 
ple precipitation  to  exceed  it  and  produce  a  negative  fluctuation  in  differential  salinity. 
In  summer  a  positive  d:fferential  salinity  pre-dominated.  Here,  the  lower  precipitation 
for  the  complex  case  allowed  the  salinity  of  this  case  to  increase. 

The  very  high  simple  precipitation  case  provided  the  best  approximation  to  ac- 
tual MLD  of  all  study  cases,  both  simple  and  complex.  The  large  amount  of  precipi- 
tation produced  consistently  shallower  MLD  s  than  the  other  simple  cases.  Additionally, 
it  provided  shallower  MLD's  than  occurred  during  the  complex  case  when  forced  by  low 
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Fig.  26.  Simple  Precipitation  Forcing  Predicted  Salinity:  (a)  dotted  line  was  low 
precipitation  forced,  (b)  dashed  line  was  medium  precipitation  (c) 
chained-dashed  was  high  precipitation  (d)  solid  line  was  very  high  precip- 
itation. 
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Fig.  27.  Subtracted  Mixed  Layer  Values  From  Higli  Complex  and  Simple 
Cases:  All  values  are  the  result  of  subtracting  complex  high  precipitation 
forcing  output  by  Low  precipitation  forcing,  (a)  differential  MLD,  (b) 
differential  temperature,  (c)  differential  salinity,  and  (d)  complex  precipi- 
tation forcinc. 
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Fig.  28.  Simple,  Complex,  and  Actual  Mixed  Layer  Values:  In  all  three  figures 
the  solid  line  represents  simple  high  forcing,  the  dotted  line  represents 
complex  forcing. 
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winter  winds  and  in  all  summer  cases  when  low  winds  occur  and  MLD's  tend  to  be 
shallowest  (Fig.  28  on  page  61). 

In  winter,  there  was  close  agreement  between  all  complex  and  simple  forced 
cases  and  actual  temperature.  In  summer,  the  shallowest  MLD's  of  the  simple  forced 
case  produced  the  highest  temperature,  which  though  best  approximated  the  actual 
temperature  in  summer,  were  still  less  than  that  of  the  actual  case  by  2.4°C. 

As  with  complex  forcing,  the  high  simple  case  provided  the  best  approximation 
to  the  actual  salinity  values.  The  simple  high  rain  case  yielded  a  larger  negative  differ- 
ential salinity  in  winter  and  a  lower  salinity  in  summer  than  the  high  complex  case,  while 
still  retaining  a  cyclic  behavior  with  a  one-year  period.  Its  cyclic  response  was  not  as 
clear  as  in  the  complex  case,  however. 
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VI.     SUMMARY 

This  study  was  conducted  to  examine  the  discrete  precipitation  event  effects  on  the 
short-term  and  seasonal  evolution  of  ocean  mixed  layer  temperature  and  salinity  struc- 
ture. Data  for  this  study  was  for  the  years  1968  and  1969  at  Ocean  Station  "P"  located 
in  the  Northeast  Pacific  Ocean.   Two  kinds  of  experiments  were  performed: 

1.  the  first  was  to  simulate  the  response  of  the  mixed  layer  to  discrete  (isolated)  pre- 
cipitation events,  and 

2.  the  second  was  a  quasi-realistic  "complex"  precipitation  forcing  experiment,  re- 
presenting a  reahVic  distribution  in  time  of  synoptic  events  over  a  13-month  pe- 
riod. 

The  hypothetical  experiment  was  conducted  using  twelve  separate  single  discrete 
precipitation  events  over  a  400-day  period.  It  examined  the  specific  effects  of  a  discrete 
precipitation  event  as  a  function  of  season  without  contamination  by  other  complex 
rainfall  history.  This  experiment  showed  that  the  effect  of  a  single  discrete  rain  event 
can  vary  with  season  and  the  MLD  (at  onset  of  the  rain  event),  as  well  as  with  the  in- 
tensity of  the  wind  stress.  Additionally,  it  showed  that  a  single  rain  event  can  have  a 
lasting  effect  on  mixed  layer  depth  and  temperature  for  up  to  55  days  after  the  event 
stops. 

For  the  main  experiment  a  "complex"  precipitation  modeling  scheme,  producing 
discrete  daily  precipitation  events  from  surface  meteorological  observations  was  devised 
and  tuned  to  yield  realistic  precipitation  input  for  the  mixed  layer  model.  The  precipi- 
tation model  was  adjusted  to  give  differing  amounts  of  total  precipitation.  Four  different 
precipitation  "intensities"  were  tested:  32  cm  400-day,  48  cm; 400-day,  96  cm,400-day, 
and  160  cm; 400-day.  The  50  cm '400-day  case  was  chosen  to  yield  the  same  net  precip- 
itation as  the  observed  precipitation  data.  These  quasi-realistic  precipitation  forcings 
were  applied  to  the  mixed  layer  model.  The  resulting  values  of  MLD,  temperature,  and 
salinity  for  this  400-day  period  were  examined  by  comparing  them  with  output  values 
when  the  model  was  forced  by  "simple"  constant  net  precipitation  with  the  same  total 
precipitation  as  in  the  complex  cases.  Finally,  the  values  predicted  by  the  model  under 
various  precipitation  forcing  conditions  were  compared  to  actual  CTD  values. 

In  all  study  cases  the  MLD's  approximated  the  MLD  cycle  of  the  actual  data.  The 
observed  MLD's  were  shallower  than  predicted  by  the  model  by  10m  or  more.  Thus,  the 
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model  runs  with  the  most  precipitation  best  approximated  the  observed  MLD's.  It  was 
noted  that  the  greatest  discrepancy  occurred  in  the  first  winter  season  of  the  400-day 
study  period.   However,  a  large  data  void  existed  during  this  time  period. 

As  with  MLD,  the  temperatures  of  these  studies  generally  approximated  the  tem- 
perature cycle  and  form  of  the  actual  data.  The  temperature  in  summer  provided  the 
greatest  discrepancy  in  temperature  magnitude  with  actual  temperatures  (exceeding 
those  of  the  output  data  by  at  least  ~  2.5°  C  in  the  complex  case  and  ~  2.5°C  in  the 
simple  case). 

Predicted  salinity  provided  the  greatest  discrepancy  between  the  modeled  and  the 
observed  cycle.  Rain  rates  of  32  cm  yr  and  48  cm/yr  (that  equal  to  the  observed  total) 
yielded  salinity  values  which  were  non-cyclic  and,  overall,  higher  in  value  than  the  ob- 
servations. The  96  cm'400-day  case  (both  complex  and  simple)  best  approximated  the 
salinity  observations.  Predicted  salinity  was  less  than  observed  in  winter  (by  ~  OAOppm 
in  the  complex  case  and  ~QA\ppm  in  the  simple  case)  and  greater  than  observed  in 
summer  (by  ~  O.Olppm  in  the  complex  case  and  ~  O.Olppm  in  the  simple  case).  The  160 
cm'400-day  precipitation  cases  produced  a  constant  freshening  of  the  mixed  layer  in 
both  winter  and  summer. 

In  the  final  analysis,  while  the  discrete  precipitation  events  had  a  different  effect  on 
the  mixed  layer  than  the  simple  cases,  the  difference  was  slight.  Increased  levels  of  pre- 
cipitation did  enhance  the  difference  between  constant  rain  and  the  quasi-realistic  rain 
cases,  as  seasonal  differences  in  precipitation  were  amplified.  The  96  cm/400-day 
quasi-realistic  precipitation  forcing  case  was  selected  as  the  most  likely  precipitation 
forcing  for  Ocean  Station  "P"  during  this  time  period.  Despite  magnitude  discrepancies 
with  observed  rain  at  Station  "P,"  the  form  and  cycle  of  the  MLD's  and  temperatures 
with  respect  to  time  were  well  approximated  by  this  level  of  precipitation.  The  most 
important  point  for  selection  of  this  precipitation  case  was,  however,  the  agreement  with 
the  actual  salinity  cycle.  Here,  the  simple  high  precipitation  case  also  approximated  the 
shape  of  the  salinity  curve,  but  was  not  as  cyclical,  especially  during  the  winter  in  1969. 
This  was  probably  because  the  complex  precipitation  case  produced  some  periods  with- 
out rain  allowing  for  shallowing,  and,  therefore,  decreased  salinity.  Even  in  winter,  while 
the  simple  cases  produced  continuous  rain.  One  such  realistically  modeled  period  acted 
to  enhance  cyclical  behavior  of  salinity  in  the  early  winter  of  1969,  adding  to  MLD 
deepening  and  adding  less  fresh  water  during  this  period  than  that  of  the  simple  case. 
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There  are  several  reasons  why  observed  mixed  layer  temperature  and  salinity  values 
may  not  agree  with  model  predicted  values: 

1.  the  inaccuracies  and  voids  in  CTD  data, 

2.  the  inaccuracies  in  model  forcing,  especially  in  producing  realistic  complex  precip- 
itation forcing, 

3.  the  inaccuracies  in  model  initial  conditions  for  T(z)  and  S(z), 

4.  the  inaccuracies  in  representing  advection,  and 

5.  the  mixed  layer  theory  inaccuracies. 

Data  voids  do  exist  in  the  available  CTD  data  at  Ocean  Station  "P."  The  largest 
data  void  occurred  during  the  first  95  days  of  the  study.  At  this  time,  there  was  an  ob- 
servation at  day  0  (Nov  23,  1968)  and  then  no  observations  until  the  95th  day  (Feb  26, 
196S).  Additionally,  data  was  sparse  from  days  95  to  180  and  days  300  to  400.  This 
would  not  be  significant  for  data  with  small  variance,  as  with  the  sea  surface  temper- 
ature. However,  for  data  with  strong  variation  such  as  the  MLD  and  salinity  data,  the 
resulting  analysis  may  be  biased  by  lack  of  information  at  times  of  critical  changes  in 
mixed  layer  evolution. 

Inaccuracies  in  prec  pitation  forcing  for  the  model  may  exist.  Four  net  precipitation 
rates,  each  multiples  of  the  first,  were  chosen  to  bracket  the  most  realistic  precipitation 
amounts  for  Station  "P."  The  fact  that  the  48  cm'400-day  case  did  not  approximate  the 
salinity  cycle  as  well  as  the  96  cm/400-day  case,  suggests  that  the  observed  precipitation 
rate  may  be  lower  than  the  actual  amount.  Allen  (1963)  outlined  the  significant  number 
of  potential  rain  error  measurement  sources  onboard  ship  at  Ocean  Station  "P."  These 
were: 

1 .  the  wind  effects, 

2.  the  shadowing  by  nearby  objects, 

3.  the  evaporation  in  the  gauge, 

4.  the  sea  spray  entering  the  gauge,  and 

5.  the  tilting  of  the  gauge  orifice  due  to  roll  and  pitch  of  the  ship. 

Additionally,  the  rain  model  scheme  used  in  this  study  was  very  crude.  It  may  have 
inaccurately  approximated  the  strength  and  duration  of  storm  events,  and  would  have 
been  completely  insensitive  to  local  convection  events  (rain  without  wind).  Such  local 
convective  events  may  account  for  a  lower  observed  salinity,  shallower  MLD's,  and 
higher  temperatures  in  summer. 
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Inaccuracies  in  meteorological  data  used  in  this  study  to  compute  heat  and  mo- 
mentum fluxes  may  also  exist.  This  data  was  collected  in  1968  and  1969,  and  collection 
techniques  have  improved  since.  In  addition  to  the  wind  stress  and  heat  flux  errors  in- 
troduced directly  into  the  model  by  this  data,  these  errors  could  produce  unrealistic 
precipitation  and  evaporation  events  as  well. 

Predicting  horizontal  advection  is  difficult,  and  given  only  surface  meteorological 
data  and  scattered  CTD  data  it  is  not  possible  with  any  accuracy.  Advective  events  have 
been  known  to  occur  at  this  station,  as  discussed  in  Chapter  I.  Such  an  event  could 
greatly  modify  mixed  layer  structure.  Advection  was  neglected  for  these  studies,  recog- 
nizing that  an  increase  in  advection  of  water  of  higher  salinity  has  the  same  effect  as 
decreasing  rain  and  is  necessary  to  balance  the  salinity  budget  in  the  longer  interannual 
term.   This  was  another  reason  for  examining  various  rain  amounts. 

It  is  recommended  that  this  experiment  be  repeated  using  observed  discrete  precipi- 
tation if  possible.  Alternatively,  Tucker  (1961)  suggested  a  rainfall  prediction  scheme 
utilizing  weather  codes,  using  the  standard  international  two-digit  system.  These 
follow-on  studies  could  help  answer  the  question  as  to  the  discrepancy  in  actual  versus 
measured  rainfall.  Finally,  precipitation  from  such  sources  as  satellite  imaging  could  be 
interpreted  as  with  Hakkarien  and  Adler  (1987)  and  applied  to  the  mixed  layer  model. 
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